AD-A201  862 


ni;  F(iE  nof' 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey ,  California 


DTIC 

SELECTEB^ 
0EC29«88|j 

q  H 

A  MICROCOMPUTER  SIMULATION  PROGRAM  TO  MODEL 
TRANSIENT  AND  STEADY-STATE  DETECTION  OF  AN 
EVADING  SUBMARINE  BY  A  SEARCHING  SUBMARINE 
IN  A  FALSE  TRANSIENT  ENVIRONMENT 

by 

J.  Brad  Kratovil 
September  1988 

Thesis  Co-Advisors:  R.  N.  Forrest 

James  N.  Eagle 


Approved  for  public  release;  distribution  Is  unlimited 


la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUM8ER(S) 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  MARKINGS 


6a  NAME  OF  PERFORMING  ORGANIZATION 

Naval  Postgraduate  School 


6b  OFFICE  SYMBOL 
(If  iPQjKiblt) 


6c.  ADDRESS  {City,  State.  *nd  /IP  Code) 

Monterey,  California  93943-5000 


8a.  NAME  OF  FUNDING  <  SPONSORING 
ORGANIZATION 


8c.  ADDRESS  (City,  Sutt,  ind  /IP  Code) 


8b  OFFICE  SYMBOL 
(If  }ppiieibi«) 


3  DISTRIBUTION  ^AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  is  unlimited 


S  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7a  NAME  OF  MONITORING  ORGANIZATION 

Naval  Postgraduate  School 


7b  ADDRESS  (Oty,  State,  a/Mt  Z/PCoete; 

Monterey,  California  93943-5000 


ft  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


to  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 


PROJECT 

NO 


WORK  UNIT 
ACCESSION  NO 


1  TITLE  (/nc/ude  Security  CFaMiFicafionj  A  Microcomputer  Simulation  Program  to  Model  Transient 
and  Steady-State  Detection  of  an  Evading  Submarine  by  a  Searching 
Submarine  in  a  Falso  Transie 


12  PERSONAL  AUTHOR(S) 


Kratovil,  j.  Brad 


13a  TYPE  OF  REPORT 

Master's  Thesis 


'3b  T  ME  COVERED 
FROM  TO 


14  DATE  OF  REPORT  (Year,  Month,  OayJ  jl5  PAGE  COUNT 

1988  September  I  128 


16  SUPPLEMENTARY  NOTATION  The  views  expressed  in  this  thesis  are  those  of  the 
author  and  do  not  reflect  the  official  policy  or  position  of  the 
Department  of  Defense  or  the  U.S.  Government _ 


18  SUBJECT  TERMS  (COQtinuP  on  rtvtnt  if  nrcemry  aripF  identify  by  (ilotk  number) 

Submarine,  Transient,  Submarine  Simulation, 
Submarine  Transient,  Submarine  Detection. 


COSATi  COOES 


SUBGROUP 


19  ABSTUifiT,  (Cononoe  on,  reverse  if  nerejuiy  end  j^ntify  J)y  block  number) 

TniS  tnesls  aescribes  the  development  and  demonstrates  the  use  of 
SUBTRAN,  a  submarine  transient  detection  computer  simulation  model. 
SUBTRAN  provides,  at  the  microcomputer  level,  a  framework  that  can  be 
used  to  Investigate  how  transient  detection  and  false  transient 
detection  opportunities  affect  the  expected  time  to  steady-state 
(continuous)  detection.  Monte  Carlo  methods  are  employed  to  simulate  a 
submarine  versus  submarine  passive  acoustic  detection  search  scenario. 
The  scenario  terminates  when  steady-state  detection  occurs.  Detection 
is  modeled  using  a  signal  excess  threshold  crossing  model.  Random 
fluctuations  in  the  acoustic  signal  excess  are  modeled  using  a  Lambda- 
"igma  Jump  process.  Both  submarines  are  assumed  to  have  a  fixed  speed 
nd  are  constrained  within  a  defined  search  area.  The  "target" 
ubmarlne  is  assumed  to  be  unable  to  counter-dete^'t  the  "searching" 
ubmarine.  Transient  signals  and  false  transient  signals  are  determined 


20  DISTRIBUTION  I  availability  OF  ABSTRACT  21  ABSTRACT  SECURITY  CLASSIFICATION 

□tuNCLASSlFlED/UNLIMlTED  □  SAME  AS  RPT  □  DTIC  USERS  Unc  la<R  S  i  f  i  pd 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL  I22b  TELEPHONE  f/nc/u<ft  AreeCode)!  c2c  OFFICE  SYMBOL 


DO  FORM  1473, 84  MAR 


83  APR  edition  may  be  used  until  exhausted 
All  other  editions  are  obsolete 


SECURITY  classification  OF  THIS  PAGE 

•  U.S.  GftvtPftftMftl  Pieties  Of*<c« 


1 


tccwmrv  euiMineATWM  mw  m«i 


19.  ABSTRACT  (cont) 


by  independent  Poisson  processes.  Sununary  statistics  for  the  times  to 
detection  are  provided  as  output  of  the  simulation. 


Copv  , 
iWspEcreo  / 

V  1 


1  Aooasslon  For  ^  | 

BTIS  GRAAI 

DTIC  TAB 

□ 

Unannounced 

□ 

Justlf loatloa. 

- — f 

- 1 

By - 

Distribution/ 

Availability  Codes 
Avail  and/or 
Dlst  Special 


SECURITY  CLAttiriCATION  Or  THIS  RACE 


Approved  for  public  release;  distribution  is  unlimited 


A  Microcomputer  Simulation  Program  to  Model 
Transient  and  Steady-State  Detection  of  an  Evading 
Submarine  by  a  Searching  Submarine  in  a  False 
Transient  Environment 


by 


J.  Brad  Kratovil 
Lieutenant,  United  States  Navy 
B.S.,  University  of  Southern  California,  1982 

Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
September  1988 


Author: 
Approved  By: 


ABSTRACT 


^  This  thesis  describes  the  development  and  demonstrates  the  use  of 
SUBTRAN,  a  submarine  transient  detection  computer  simulation  model. 
SUBTRAN  provides,  at  the  microcomputer  level,  a  framework  that  can  be  used 
to  investigate  how  transient  detection  and  false  transient  detection  opportunities 
affect  the  expected  time  to  steady-state  (continuous)  detection.  Monte  Carlo 
methods  are  employed  to  simulate  a  submarine  versus  submarine  passive 
acoustic  detection  search  scenario.  The  scenario  terminates  when  steady-state 
detection  occurs.  Detection  is  modeled  using  a  signal  excess  threshold  crossing 
model.  Random  fluctuations  in  the  acoustic  signal  excess  are  modeled  using  a 
Lambda-Sigma  Jump  process.  Both  submarines  are  assumed  to  have  a  fixed 


speed  and  are  constrained  within  a  defined  search  area.  The  target”  submarine 
is  assumed  to  be  unable  to  counter-detect  the  searching  submarine.  Transient 
signals  and  false  transient  signals  are  determined  by  independent  Poisson 
processes.  Summary  statistics  for  the  times  to  detection  are  provided  as  output 


of  the  simulation.  ^ 
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I.  INTRODUCTION 


Passive  detection  of  a  submarine  can  be  divided  into  two  mutually  exclusive 
events:  steady-state  detection  and  transient  detection.  Steady-state  detection 
involves  the  passive  acoustic  monitoring  of  a  submarine’s  continuous  operating 
sounds.  With  today’s  submarine  quieting  capabilities,  steady-state  detection 
typically  requires  short  ranges.  On  Uie  other  hand,  transient  detection  is  the 
detection  of  a  submarine’s  intermittent  sounds  that  are  produced  during 
infrequent  operations.  Due  to  the  high  signal  levels  generated  by  many  of  these 
infrequent  operations,  transient  detection  can  occur  at  very  long  ranges. 

Of  great  concern  to  an  evading  submarine  is  the  ability  to  freely  move  about 
without  its  actions  being  monitored  by  a  searching  submarine.  This  monitoring  is 
generally  based  on  an  evading  submarine’s  continuous  operating  sounds. 

Due  to  the  typical  short  ranges  involved,  the  chance  of  a  searching 
submarine  being  able  to  monitor  an  evading  submarine’s  continuously  generated 
signals  is  quite  small  when  a  large  search  area  is  involved.  However,  an  evading 
submarine’s  transient  signals  offer  long-range  detection  opportunities.  The 
exploitation  of  these  opportunities  may  improve  the  ability  of  the  searching 
submarine  to  gain  steady-state  detection. 

A  drawback  to  the  use  of  long  range  detection  opportunities  occurs  in  a  false 
transient  environment.  In  any  search  scenau-io,  false  detections  can  lead  the 
searcher  away  from  the  actual  target  position. 

The  purpose  of  this  thesis  is  to  develop  a  means  with  which  to  address  the 
issue  of  whether  or  not  transient  detection  opportunities  can  be  utilized  by  a 
searching  submarine  in  a  false  transient  environment  to  increase  its  chance  of 
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gaining  steady-state  detection  on  a  target  submarine.  This  will  be  accomplished 
by  use  of  a  Monte  Carlo  simulation,  SUBTRAN. 
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U.  BACKGROUND 


A.  GENERAL  DESCRIPTION  OF  SIMULATION  MODEL 

SUBTRAN  is  a  discrete  event  simulation  of  two  submarine  platforms  in  a 
area  search  scenario.  One  platform  is  designated  the  "target"  and  one  the 
"searcher".  During  each  discrete  event  time  step,  the  time  to  the  next  discrete 
event  is  computed,  steady-state  detection  is  assessed,  target  and  searcher 
positions  are  updated  and  the  simulation  time  is  advanced.  A  simulation  run  is 
completed  when  steady-state  detection  occurs  or  the  user-defined  maximum 
search  time  is  reached.  Propagation  loss  curves  and  figure  of  merits  are  used  to 
derive  the  deterministic  element  of  signal  excess.  A  Lambda-Sigma  Jump  process 
is  used  to  describe  stochastic  variations  in  signal  excess,  and  a  threshold 
crossing  model  is  used  to  evaluate  total  signal  excess  for  the  occurrence  of 
detection.  Target  characteristics  incorporated  in  the  simulation  model  include: 
target’s  course,  target’s  speed,  and  target’s  transient  occurrence  rate.  Searcher 
characteristics  incorporated  in  the  simulation  model  include  the  searcher’s  course 
emd  speed.  Target  depth,  relative  to  the  searcher,  is  assumed  constant.  Fedse 
transients  are  modeled  in  the  same  fashion  as  target  transients  by  using  a 
transient  occurrence  rate. 

B.  SEARCH  SCENARIO 

The  search  scenario  modeled  in  SUBTRAN  is  a  parallel  sweep  search  of  a 
square  area  whose  boundaries  are  known  by  both  the  target  and  searcher.  At 
the  beginning  of  each  simulation  replication,  the  target  is  positioned  such  that 
its  X  component  of  position  and  Y  component  of  position  are  independent 


.a 


uniform  random  variables  on  the  interval  (O.L).  L  is  the  length  dimension  of  the 
search  area  in  uautical  miles.  The  searcher’s  position,  at  the  beginning  of  each 
replication,  is  in  the  lower  left-hand  corner  of  the  search  area.  In  order  to 
maximize  steady-state  search  effort,  the  searcher  is  positioned  a  distance  equal 
to  one-half  a  track  width  from  both  the  bottom  and  left-hand  search  area 
boundaries.  The  search  area  is  oriented  such  that  the  top  area  boundary  is  to 
the  north  for  course  reference.  A  typical  search  area,  annotated  with  initial 
searcher  and  target  positions,  is  displayed  in  Figure  2.1. 


Nt 

X 

TARGET'S  POSITION 

SEARCHER'S  POSITION 

TM/t 

X. 

TH/2 

Figure  2.1.  T3^ical  Search  Area. 
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C.  DETECTION  MODEL 

The  basic  detection  model  used  in  SUBTRAN  converts  user-defined 
propagation  loss  curves,  figure-of-merits  (FOM),  and  Lambda-Sigma  Jump 
process  parameters  into  an  assessment  of  detection  occurrence  by  evaluating 
total  signal  excess  using  a  threshold  crossing  detection  model.  Two  tsrpes  of 
detection  are  modeled:  steady-state  detection  and  transient  detection.  Due  to  the 
current  trend  in  submarine  quieting,  convergence  zone  detection  of  the  target’s 
steady-state  sound  signature  is  assumed  not  to  be  possible.  Convergence  zone 
detection  of  the  target’s  transient  noise  signature  is  assumed  to  be  possible  but 
only  to  the  third  convergence  zone. 

1.  Propagation  Loss 

The  propagation  loss  experienced  by  sound  as  it  travels  through  the 
water  is  determined  by  many  factors.  Frequency  of  the  sound,  receiver  depth, 
source  depth,  ocean  bottom  tjrpe,  and  sound  velocity  profile  all  have  an  effect  on 
the  intensity  of  the  sound  received  from  an  underwater  source.  In  the  area 
search  scenario  modeled  in  SUBTRAN,  most  of  these  factors  remain  relatively 
static  with  the  exception  of  the  frequency  of  the  sound.  Since  steady-state  and 
transient  sounds  are  produced  differently,  they  are  likely  to  have  different 
frequency  components.  As  a  result,  two  propagation  loss  curves  are  modeled  in 
SUBTRAN,  a  steady-state  propagation  loss  curve  and  a  transient  propagation 
loss  curve. 

SUBTRAN  allows  the  user  to  input  up  to  50  values  of  range  and 
corresponding  propagation  loss  for  each  of  the  two  types  of  sound.  From  this 
input,  the  model  constructs  propagation  loss  curves  using  straight  line 
approximation  between  successive  points.  One  curve  is  constructed  for  steady- 
state  signal  loss  and  one  is  constructed  for  transient  signal  loss.  Since 


convergence  zone  detection  of  the  target’s  steady-state  sound  signature  is 
assumed  not  to  be  possible,  only  the  direct  path  remge  on  the  steady-state 
propagation  loss  curve  is  used. 

2.  Figure-of'Merit 

FOM,  defined  for  passive  sonar  as  the  maximum  propagation  loss  which 
still  results  in  a  50  percent  probability  of  signal  detection  [Ref.  l:p.  49],  is 
modeled  in  SUBTRAN  tinder  the  following  assumptions  and  conditions: 

-  The  target  cannot  counter-detect  the  searcher. 

-  The  speeds  of  the  target  and  searcher  are  held  constant. 

-  There  is  only  one  type  of  steady-state  sound  source  emitted  by  the  target 
and  it  is  omni-directional. 

-  There  is  only  one  type  of  transient  sound  source  emitted  by  the  target  and 
it  is  omni-direction^. 

•  The  ambient  noise  level  is  constant  throughout  the  simulation. 

Under  these  assumptions  and  conditions,  target  FOM  can  be 
summarized  by  two  constant  values:  steady-state  FOM  and  transient  FOM. 
SUBTRAN  requires  the  user  to  determine  the  values  of  the  steady-state  and 
transient  FOM  for  the  target  based  on  the  search  scenario  being  simulated. 

3.  Signal  Excess 

Signal  excess  (SE)  is  modeled  in  SUBTRAN  as  consisting  of  a 
deterministic  part  and  a  stochastic  part.  The  deterministic  part  of  SE  is  defined 
as  the  algebrcuc  difference  between  FOM  and  propagation  loss.  The  stochastic 
part  of  signal  excess  is  used  to  describe  the  random  fluctuations  that  occur  in 
signal  excess  with  respect  to  time.  A  Lambda-Sigma  Jump  process  is  used  to 
model  these  random  fluctuations. 
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4.  Lambda>Sigma  Jump  Process 

A  Lambda-Sigma  Jump  process  [Ref.  2:p.  8]  is  used  to  describe  the 
random  fluctuations  in  signal  excess.  This  process  models  the  time  between 
fluctuations  as  well  as  the  magnitude  of  the  fluctuations.  For  the  Lambda-Signal 
Jump  process,  the  time  between  fluctuations  is  exponentially  distributed  with 
rate  parameter  lambda  and  the  magnitude  of  the  fluctuations  is  normally 
distributed  with  mean  zero  and  standard  deviation  sigma. 

The  user  inputs  values  for  the  parameters  lambda  and  sigma. 
Appropriate  values  for  these  parameters  are  subjective  and  must  be  based  on  the 
experience  of  the  user.  A  guide  in  the  estimation  of  lambda  and  sigma,  based  on 
various  sources,  is  [Ref  3  and  Ref.  4]: 

-  Lambda  varies  between  2  and  60  hours  ’. 

-  Sigma  varies  between  3  and  9  dB. 

Other  sources  containing  information  helpful  in  determining  appropriate  values 
for  lambda  and  sigma  are  listed  by  Tehan  [Ref  4]. 

5.  Threshold  Crossing  Model 

A  threshold  crossing  model  [Ref  2:p.  7]  is  the  final  step  in  the  overall 
effort  to  model  detection.  This  model  predicts  that  detection  occurs  when  signal 
excess  is  equal  to  or  exceeds  some  given  threshold  level.  The  most  common 
threshold  level,  and  the  one  employed  in  SUBTRAN,  is  zero.  In  other  words, 
detection  occurs  whenever  signal  excess  is  greater  than  or  equal  to  zero. 

6.  Detection  Model  Employment 

The  detection  model  outlined  in  the  previous  sections  is  implemented  in 
SUBTRAN  by  extracting  the  current  detection  range  from  the  propagation  loss 
curve  and  comparing  it  to  the  range  between  the  target  and  the  searcher.  If  the 


detection  range  is  greater  than  or  equal  to  the  range  between  the  target  and  the 
searcher,  detection  occurs. 

The  current  detection  range  is  determined  by  considering  an  adjusted 
value  of  the  FOM.  This  adjusted  value  is  obtained  by  adding  the  user-defined 
value  of  FOM  and  the  current  signal  excess  fluctuation  value.  The  range  at 
which  this  adjusted  FOM  value  and  the  propagation  loss  curve  are  equal  is  the 
current  detection  range. 

7.  Steady-State  Detection  Aasesament 

SUBTRAN  is  a  discrete  event  simulation.  A  discrete  event  is  defined  as 
any  event  whose  occurrence  could  change  the  ability  of  the  searcher  to  gain 
steady-state  detection  on  the  target.  The  discrete  events  modeled  in  SUBTRAN 
are:  target  course  change,  target  boundary  reflection,  searcher  course  change, 
target  transient  signal  occurrence,  false  transient  signal  occurrence,  and  signal 
excess  fluctuation  occurrence. 

Since  SUBTRAN  is  a  discrete  event  simulation,  steady-state  detection 
must  be  assessed  for  the  time  period  until  the  next  discrete  event.  To  do  this, 
the  current  values  of  searcher  and  target  course  and  position  are  used  together 
with  the  user-defined  speed  values  to  define  the  positions  of  the  searcher  and 
target  at  any  future  time.  These  positions  are  then  compared  to  find  the  time  it 
would  take  the  range  between  the  searcher  and  target  to  decrease  to  a  value 
less  than  or  equal  to  the  current  steady-state  detection  range.  If  the  time 
calculated  is  less  than  or  equal  to  the  time  to  the  next  discrete  event,  steady- 
state  detection  occurs.  Appendix  E  contains  the  steady-state  detection  assessment 
calculations. 
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8.  Transient  Detection  Assessment 

Detection  of  target  transients  by  the  searcher  is  assessed  in  a  different 
manner  than  steady-state  detection.  The  target’s  transient  sound  signature  is 
intermittent  and  when  produced,  is  assumed  to  occtu*  instantaneously.  As  a 
result  of  these  characteristics,  transient  detection  is  assessed  in  SUBTRAN  by 
comparing  searcher  and  target  positions  at  the  time  of  a  transient  occurrence.  If 
the  range  is  within  the  transient  detection  range,  transient  detection  occurs. 

Transient  detection  range  is  extracted  from  the  transient  propagation 
loss  curve  using  the  transient  FOM.  Since  transient  detection  can  occur  out  to 
the  third  convergence  zone,  the  transient  detection  range  consists  of  four  range 
brackets:  a  direct  path  range  bracket  and  three  convergence  zone  annulus  range 
brackets. 

D.  TARGET  CHARACTERISTICS 

The  target  characteristics  modeled  in  SUBTRAN  include:  speed,  course,  and 
transient  emission  rate.  Target  depth  relative  to  the  searcher  is  assumed  to  be 
constant.  In  order  to  model  its  worst  case  vulnerability  to  detection,  the  target 
does  not  have  the  ability  to  detect  and  thus  avoid  the  searcher. 

1.  Target  Speed 

Target  speed  is  held  constant  throughout  the  simulation  at  a  value 
specified  by  the  user. 

2.  Target  Course 

The  search  scenario  modeled  in  SUBTRAN  involves  a  randomly  moving 
target  which  is  constrained  to  a  defined  search  area  by  boundary  reflection.  To 
incorporate  differing  target  course  and  time  to  course  change  distributions, 


SUBTRAN  uses  the  inverse  probability  integral  transformation  [Ref.  5:pp.  20- 
21a}  to  model  target  movement. 


The  inverse  probability  integral  transformation  method  involves 
selecting  values  of  a  desired  random  variable  whose  cumulative  distribution 
function  (CDF)  is  known  by  using  random  variables  drawn  from  a  uniform  (0,1) 
distribution.  This  method  is  illustrated  in  Figure  2.2.  A  random  variable  drawn 
from  a  uniform  (0,1)  distribution  is  entered  on  the  Y-azis  of  the  CDF 
representing  the  random  variable  desired.  The  corresponding  X-axis  value  is  the 
associated  random  variable  of  the  given  distribution. 

SUBTRAN  allows  the  user  to  input  up  to  50  values  for  both  the  target’s 
course  CDF  and  the  time  to  course  change  CDF.  The  corresponding  CDF  curves 
are  constructed  using  straight-line  approximation  between  successive  points. 

Due  to  the  randomness  of  the  target’s  movement,  a  procedure  exists  to 
constrain  the  target  to  the  search  area.  Using  its  current  course  and  speed, 
the  position  of  the  target  is  evaluated  for  the  period  until  its  next  course  change. 
If  dxuring  this  time  the  target  encounters  the  search  area  boundary,  its  course  is 
reflected  back  into  the  search  area.  The  angle  of  reflection  is  equed  to  the  angle 
of  incidence  between  the  search  area  boundary  and  the  target’s  incoming  course. 
An  example  of  target  boundary  reflection  is  shown  in  Figure  2.3. 

3c  Target  Transients 

Tsu-get  transient  occurrence  is  modeled  in  SUBTRAN  as  a  Poisson 
process  with  a  user-defined  rate  parameter.  Only  one  type  of  transient  sound  is 
assumed  to  be  produced  by  the  target.  In  addition,  when  the  transient  does 
occur,  it  is  assumed  to  be  instantaneous. 
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Figure  2.2  Inverse  Probability  Integral  Transformation  Method. 
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E.  SEARCHER  CHARACTERISTICS 

Searcher  characteristics  modeled  in  SUBTRAN  include;  speed  and  course. 
The  depth  of  the  searcher  relative  to  the  target  is  assumed  constant. 

1.  Searcher  Speed 

Searcher  speed  is  held  constant  throughout  the  simulation  at  a  value 
specified  by  the  user. 

2.  Searcher  Course 

Searcher  movement,  unlike  target  movement  is  deterministic.  The 
movement  modeled  in  SUBTRAN  consists  of  a  general  search  movement  and  a 
transient  prosecution  movement. 

The  general  search  movement  is  an  area  search  using  parallel  sweeps. 
The  spacing  between  sweeps  is  defined  as  the  track  width,  TW.  TW  can  be  user 
specified  or  defaulted.  When  defaulted,  TW  is  set  equal  to  twice  the  value  of  the 
mean  steady-state  detection  range.  To  maximize  the  amount  of  search  effort 
applied  to  the  search  area,  the  searcher  is  programmed  to  stay  a  distance  of 
(0.5)  X  (TW)  from  the  area  boundary.  This  distance  constraint  forces  an 
adjustment  to  the  searcher’s  movement  as  it  approaches  the  end  of  the  area 
search.  Specifically,  after  each  vertical  leg  of  the  search,  the  position  of  the 
search  is  checked.  If  the  horizontal  distance  between  the  searcher  and  the  search 
area  boundary  that  lies  ahead  of  the  searcher  is  less  than  (1.5)  x  (TW),  the  next 
horizontal  leg  of  the  search  is  adjusted.  The  adjustment  that  occurs  is  a 
shortening  of  the  distance  the  searcher  normally  travels.  The  normal  distance  of 
the  horizontal  leg  of  the  search  is  TW.  The  adjusted  distance,  TW^j,  is  given  by 
the  following  expression: 

TW^j  =  TW  X  INT(I/rW) 
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where, 


L  =  length  dimension  of  the  search  area 

INTdvTW)  =  integer  part  of  L/TW 

Figure  2.4  illustrates  the  track  width  adjustment  and  displa}^  the 
resulting  return  area  search  that  occurs.  Maximum  search  area  coverage  is 
achieved  when  l/TW  is  an  integer  value. 

The  prosecution  phase  of  searcher  movement  involves  suspending  or 
terminating  the  current  search,  investigating  a  transient  detection,  and  returning 
to  the  general  area  search  if  an  additional  transient  detection  does  not  occur. 

Upon  detection  of  a  transient,  the  searcher  either  suspends  the  general 
area  search  or  terminates  a  prosecution  movement,  depending  on  what  search 
phase  is  being  conducted,  and  commences  a  new  prosecution  movement.  A 
prosecution  movement  involves  steering  the  course  on  which  the  transient  was 
detected  until  the  search  area  boundary  is  approached.  When  the  distance 
between  the  searcher  and  the  search  area  boundary  along  the  searcher’s  track 
reduces  to  (0.5)  x  (TW),  the  searcher  reverses  direction  and  continues  down  the 
reciprocal  bearing  until  he  returns  to  the  point  from  which  the  prosecution 
movement  was  commenced.  From  this  point,  the  searcher  returns  to  the  position 
where  the  general  area  search  was  suspended  and  resumes  the  general  area 
search.  If  at  any  time  during  the  prosecution  movement  phase  another  transient 
is  detected,  the  searcher  terminates  the  current  movement  phase  and  commences 
a  new  one  based  on  the  most  recent  transient  detection. 


13 


Figure  2.4.  Searcher  Track  Width  Adjustment 
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F.  FALSE  TRANSIENT  MODEL 

False  transients  are  included  in  SUBTRAN  to  represent  transient  sound 
sources  that  the  searcher  may  mistakenly  identify  as  target  transients.  The 
characteristics  of  false  transients  that  are  modeled  in  SUBTRAN  include:  the 
rate  of  occurrence,  the  position  at  time  of  occurrence,  and  the  ability  of  the 
searcher  to  detect  the  occurrence. 

False  transient  occurrence  is  modeled  as  a  Poisson  process.  As  a  result,  the 
rate  of  false  transient  occurrence  is  exponentially  distributed  with  a  user-dehned 
rate  parameter. 

In  an  effort  to  best  represent  natural  conditions,  the  position  from  which  a 
false  transient  is  emitted  is  assumed  to  be  uniformly  distributed  over  the  search 
area. 

The  ability  of  the  searcher  to  detect  the  occurrence  of  a  false  transient  is 
assessed  in  the  same  manner  as  the  detection  of  real  transients.  SpeciBcally,  the 
range  between  the  searcher  and  the  position  that  the  false  transient  was  emitted 
from  is  compared  to  the  detection  ranges  extracted  from  the  transient 
propagation  loss  curve.  The  FOM  for  false  transients  is  assumed  to  be  the  same 
as  the  FOM  for  real  transients. 

G.  SUMMARY  OF  IMPORTANT  MODEL  ASSUMPTIONS 

The  migor  assumptions  made  in  the  development  of  SUBTRAN  include: 

-  The  search  area  is  square  and  its  boundaries  are  known  to  both  the  target 

and  searcher. 

-  No  convergence  zone  steady-state  detection. 

-  Third  convergence  zone  transient  detection. 

-  Target  cannot  counter-detect  the  searcher. 

-  Only  one  t3T)e  of  target  transient  sound. 
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Target  transient  sounds  are  instantaneous. 

FOM  is  constant. 

Random  fluctuations  in  signal  excess  are  described  by  a  Lambda-Sigma 
Jump  process  having  the  values  of  lambda  and  sigma  specified  by  the  user. 

Target  and  false  transients  are  independent  Poisson  processes. 

False  transients  sounds  have  the  same  FOM  as  target  transient  sounds. 

Target  position  is  constrained  to  the  search  area  by  reflection. 
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HL  DESCRIPnON  OF  SUBTRAN 


A-  PROGRAM  CODE  AND  STRUCTURE 

SUBTRAN  is  coded  in  Fortran  77  [Ref.  6]  and  compiled  using  IBM 
Professional  Fortran.  It  is  designed  to  be  implemented  on  a  personal  computer. 
The  program  is  structured  in  modular  form  with  major  tasks  such  as  boundary 
reflection,  assessment  of  steady-state  detection,  and  searcher  movement  being 
accomplished  by  separate  subroutines.  This  t3rpe  of  structuring  facilitates  easy 
alteration  of  major  program  functions  by  a  knowledgeable  programmer.  However, 
a  compiled  version  of  the  program  is  all  that  is  needed  by  most  users.  The 
compiled  version  was  designed  to  allow  a  user,  who  accepts  the  assumptions 
made  in  this  thesis,  to  change  a  select  set  of  input  parameters  in  an  input  data 
file  and  not  have  to  perform  any  programming. 

B.  MAIN  PROGRAM  AND  ASSOCIATED  SUBROUTINES 

The  programming  code  for  SUBTRAN  is  broken  up  into  19  parts:  a  main 
program  and  18  supporting  subroutines.  The  interrelationship  of  the  various 
parts  is  shown  in  Figure  3.1. 

1.  Main  Program 

The  main  program,  SUBTRAN. FOR,  controls  the  call  of  the  various 
subroutines.  Specific  functions  performed  in  SUBTRAN.FOR  include: 
initiaUzation  of  counters,  initial  positioning  of  the  searcher  and  target, 
calculation  of  the  time  to  the  next  discrete  event,  control  of  the  simulation  time 
clock,  storage  of  the  times  to  steady-state  detection,  and  control  of  the  number  of 
simulation  replications. 
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2.  Input  Subroutine 

The  input  subroutine,  IN.FOR,  is  responsible  for  reading  the  input  data 
file  and  passing  the  input  parameters  to  the  main  program.  Other  functions 
performed  in  IN.FOR  include:  printing  the  input  parameters  to  the  user’s 
terminal  for  verification  and  printing  the  input  parameters  to  the  output  data 
file  SUBTRAN.OUT. 

3.  Propagation  Loss  Subroutine 

The  propagation  loss  subroutine,  PLOSS.FOR,  is  responsible  for 
calculating  the  steady-state  and  transient  detection  ranges.  Major  parameters 
passed  to  PLOSS.FOR  by  the  main  program  include:  the  propagation  loss  curve 
values  and  the  adjusted  FOM. 

4.  Target  Course  Subroutine 

The  target  course  subroutine,  TGCR.FOR,  is  responsible  for  calculating 
the  target’s  next  course  and  time  to  course  change.  Major  parameters  passed  to 
TGCR.FOR  by  the  main  program  include:  the  target’s  course  and  time  to  course 
change  CDF  curve  values  and  a  random  number  seed. 

5.  Target  Transient  Subroutine 

The  target  transient  subroutine,  TGTR.FOR,  is  responsible  for 

calculating  the  time  to  the  target’s  next  transient.  Major  parameters  passed  to 
TGTR.FOR  by  the  main  program  include:  the  target’s  transient  occurrence  rate 
parameter  and  a  random  number  seed. 

6.  Searcher  Course  Subroutine 

The  searcher  course  subroutine,  SRCR.FOR,  is  responsible  for 

calculating  the  searcher’s  next  course  and  time  to  course  change  during  general 
area  search  movement.  Miyor  parameters  passed  to  SRCR.FOR  by  the  main 
program  include:  the  searcher’s  position  and  speed,  track  width,  length 
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dimension  of  the  search  area,  and  searcher’s  course  flag.  The  searcher’s  course 
flag  indicates  which  part  of  the  general  area  search  is  to  be  performed  next. 

7.  Transient  Course  Subroutine 

The  transient  course  subroutine,  TRCR.FOR,  is  responsible  for 
calculating  the  searcher’s  next  cotirse  and  time  to  course  change  during  transient 
prosecution  movement.  Major  parameters  passed  to  TRCR.FOR  by  the  main 
program  include:  the  searcher’s  position  and  speed,  track  width,  length 
dimension  of  the  search  area,  and  searcher’s  course  flag  and  transient  course 
flag.  The  searcher’s  transient  course  flag  indicates  two  things:  whether  the 
searcher  is  performing  a  general  area  search  movement  or  transient  prosecution 
movement,  and  if  a  transient  prosecution  movement  is  being  performed,  which 
leg  of  the  movement  is  to  be  performed  next. 

8.  Steady*State  Detection  Subroutine 

The  steady-state  detection  subroutine,  SSDET.FOR,  is  responsible  for 
calculating  the  time  to  steady-state  detection.  Miyor  parameters  passed  to 
SSDET.FOR  by  the  main  program  include:  the  position  of  the  searcher  and 
target,  the  course  and  speed  of  the  searcher  and  target,  and  the  steady-state 
detection  range. 

9.  ’Transient  Detection  Subroutine 

The  transient  detection  subroutine,  TRDET.FOR,  is  responsible  for 
assessing  the  occurrence  of  target  transient  and  false  transient  detection.  Migor 
parameters  passed  to  TRDET.FOR  by  the  main  program  include:  the  position  of 
the  searcher  and  the  transient  emission,  and  the  transient  detection  ranges. 

10.  Lambda-Sigma  Jump  Subroutine 

The  Lambda-Sigma  Jump  subroutine,  LSJ.FOR,  is  responsible  for 
calculating  the  magnitude  of  the  next  signal  excess  fluctuation  as  well  as  the 
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time  to  occurrence.  Major  parameters  passed  to  LSJ.FOR  by  the  main  program 
include:  the  Lambda-Sigma  Jump  error  process  rate  parameter  and  standard 
deviation,  and  a  random  number  seed. 

11.  Falae  Transient  Subroutine 

The  false  transient  subroutine,  FALTR.FOR,  is  responsible  for 
calculating  the  position  of  and  time  to  the  next  false  transient.  Major  parameters 
passed  to  FALTR.FOR  by  the  main  program  include;  the  false  transient  rate 
parameter,  the  length  dimension  of  the  search  area,  and  a  random  number  seed. 

12.  Position  Subroutine 

The  position  subroutine,  POSIT.FOR,  is  responsible  for  updating  the  X 
and  Y  positions  of  the  searcher  and  target.  Major  parameters  passed  to 
POSIT.FOR  by  the  main  program  include:  the  previous  X  and  Y  positions  of  the 
searcher  and  target,  the  course  and  speed  of  the  searcher  and  target,  and  the 
time  passed  sioce  the  last  update. 

13.  Boundary  Reflection  Subroutine 

The  boundary  reflection  subroutine,  REFLCT.FOR,  is  responsible  for 
calculating  the  targets  next  boundary  reflection  course  and  the  time  to  reflection. 
Major  parameters  passed  to  REFLCT.FOR  by  the  main  program  include;  the 
target’s  position,  the  target’s  course  and  speed,  and  the  length  dimension  of  the 
search  area.  REFLECT.FOR  is  used  by  the  TRACK.FOR  subroutine  in  addition 
to  the  main  program. 

14.  Track  Subroutine 

The  track  subroutine,  TRACK.FOR,  is  responsible  for  calculating  the 
course  and  time  necessary  for  the  searcher  to  reach  a  given  point  in  the  search 
area.  Major  parameters  passed  to  TRACK.FOR  by  the  main  program  include:  the 
positions  of  the  searcher  and  the  point  to  be  reached,  the  speed  of  the  searcher, 


and  the  length  dimension  of  the  search  area.  TRACK.FOR  is  used  by  the 
TRCR.FOR  subroutine  in  addition  to  the  main  program. 

15.  Output  Subroutine 

The  output  subroutine,  OUT.FOR,  is  responsible  for  writing  the  output 
data  to  the  user’s  terminal  and  to  the  two  output  data  files. 

16.  Miscellaneous  Subroutines 

Four  miscellaneous  subroutines  are  called  by  the  meun  program  and 
other  subroutines:  HSCRN.FOR,  LRNDPC.FOR,  LNORPC.FOR,  and 
LGAMPC.FOR. 

The  subroutine  HSCRN.FOR  is  used  by  the  main  program,  the  input 
subroutine,  and  the  output  subroutine  to  stop  the  output  display  on  the  user’s 
terminal. 

The  subroutines  LRNDPC.FOR,  LNORPC.FOR,  and  LGAMPC.FOR  are 
used  by  the  main  program  and  various  subroutines  to  generate  random  numbers. 
All  three  are  an  adapted  version  of  the  subroutines,  with  the  same  names, 
developed  by  Lewis,  Orav,  and  Uribe  [Ref.  7]. 

LRNDPC.FOR  is  used  to  generate  uniform  (0,1)  rtmdom  variables, 
LNORPC.FOR  is  used  to  generate  normtd  (0,1)  random  variables,  and 
LGAMPC.FOR  is  used  to  generate  exponential  (1)  random  variables. 

C.  PROGRAM  OUTPUT 

The  output  provided  by  SUBTRAN  consists  of:  a  verification  of  input 
parameters,  various  output  statistics,  and  the  specific  times  to  detection.  The 
output  is  written  to  two  separate  output  files  named  SUBTRAN.OUT  and 
TIME.OUT. 


1.  Verification  of  Input  Parameters 

All  parameters  entered  via  the  input  data  file  are  printed  in  the  output 
file  SUBTRAN.OUT.  This  allows  user  verification  of  input  as  well  as  a  hardcopy 
display. 


2.  Various  Output  Statistics 

SUBTRAN  performs  elementary  statistical  analysis  on  the  times  to 
detection  and  provides  the  results  in  the  output  data  file  SUBTRAN.OUT.  The 
statistics  include:  mean  and  standard  deviation  of  the  times  to  detection,  95% 
confidence  bound  for  the  mean,  estimate  of  the  mean  time  to  detection  assuming 
the  times  to  detection  are  exponentially  distributed,  R^  confidence  factor  for  the 
estimated  mean  time  to  detection,  and  the  cumulative  distribution  function  of 
the  times  to  detection  in  tabular  form.  Calctilations  are  performed  on  the  times 
to  detection  conditioned  on  no  detection  at  time  zero.  In  other  words,  runs  with 
detection  time  zero  were  disregarded.  The  number  of  simulation  runs  with  a  zero 
detection  time  is  included  in  the  output  summary. 

The  mean  and  standard  deviation  are  computed  using  normal  methods. 
For  the  replications  in  which  steady-state  detection  does  not  occur,  T„„  is  used 
as  the  detection  time  for  the  calculation  of  the  mean  and  standard  deviation. 
This  results  in  a  lower  bound  on  the  actued  mean  time  to  detection  and  standard 
deviation.  The  number  of  times  to  detection  that  are  set  to  T„„  is  provided  in 
the  output  summary. 

The  95%  confidence  bound  for  the  mean  is  calculated  using: 

95%  bound  =  o/N" 


where. 


o  =  standard  deviation  of  the  times  to  detection. 
N  =  number  of  simulation  replications. 


] 

j 

I 


I 


I 


23 


The  estimated  mean  time  to  detection,  under  the  assumption  that  the 
times  to  detection  are  exponentially  distributed,  is  calculated  using  straight  line 
regression  of  the  natural  log  of  the  empirical  survival  function  on  the  ordered 
times  to  detection  [Ref.  8:pp.  604-607].  The  stirvival  function  of  an  exponential 
distribution  when  plotted  on  a  natural  log  scale  is  a  straight  line.  The 
magnitude  of  the  slope  of  the  straight  line  is  equal  to  the  rate  parameter  of  the 
exponential  distribution.  The  reciprocal  of  the  rate  parameter  is  the  estimated 
mean  time  to  detection.  Replications  in  which  steady-state  detection  did  not 
occwr  are  disregarded  in  the  regression  under  the  assumption  that  if  the 
distribution  is  truly  exponential,  their  inclusion  would  not  affect  the  estimate  of 
the  slope. 

The  R*  confidence  factor  for  the  estimated  mean  time  to  detection 
describes  the  proportion  of  the  variation  in  the  sample  values  of  the  natural  log 
of  the  survival  function  that  the  regression  estimate  accounts  for.  R*  was 
calculated  using  normal  methods  [Ref.  8:p.  644]. 

The  cumulative  distribution  function  of  the  times  to  detection  is 
empirically  derived  from  the  ordered  times  to  detection  [Ref.  9:pp.  11-16]  and  is 
provided  in  tabular  form.  An  example  of  the  output  statistics  provided  by 
SUBTRAN  in  contained  in  Appendix  D. 

3.  Specific  Times  to  Detection 

The  specific  times  to  detection  are  provided  in  the  output  data  file 
named  TIME.OUT.  The  times  are  ordered  from  shortest  to  longest  with  zero 
times  removed  and  T„„  value  included  as  necessary. 
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IV.  VALIDATION  OF  SUBTRAN 


Validation  of  a  simulation  model  can  be  broken  up  into  two  major  parts: 
validation  that  the  computer  code  is  functioning  as  expected  and  validation  that 
the  simulation  model  is  an  accurate  assessment  of  the  real  world  situation  being 
modeled. 

With  regard  to  SUBTRAN,  validation  that  the  computer  code  is  functioning 
as  expected  was  conducted  using  two  separate  techniques.  First,  the  simulation 
was  conducted  with  a  fixed  target  and  the  time  to  detection  was  compared  to  a 
uniform  distribution.  Second,  the  simulation  was  conducted  with  a  fixed  searcher 
and  the  time  to  detection  was  compared  to  model  results  for  a  random  tour 
target  developed  by  Abd  El-Fadeel  [Ref.  10].  The  results  of  these  two  tests  were 
used  to  validate  the  proper  simulation  coding  of  searcher  movement,  target 
positioning,  target  movement,  target  boimdary  reflection,  and  steady-state 
detection  assessment. 

Areas  of  program  code  not  rigorously  validated  include  searcher  transient 
prosecution,  transient  occurrence,  and  false  transient  occurrence.  In  the  case  of 
transient  occurrence  and  false  transient  occurrence,  the  program  code  is 
relatively  straightforward  and  rigorous  validation  is  not  necessary.  The  progreun 
code  for  searcher  transient  prosecution  is  more  complex  and  validation  was 
conducted  by  reviewing  searcher  position  data  after  simulation  runs  in  which 
transients  were  allowed.  The  data  reviewed  indicated  the  searcher  was 
performing  transient  prosecutions  as  modeled. 

Validation,  in  a  strict  sense,  that  the  simulation  model  is  an  accurate 
assessment  of  the  real  world  situation  being  modeled  is  not  possible.  There 
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exists  little  to  no  data  that  cotild  be  used  in  a  rigorous  validation  of  SUBTRAN. 
However,  the  results  produced  by  the  ^cample  runs  of  SUBTRAN  lend  credibility 
to  the  modeling  technique  used.  The  results  clearly  indicated  intuitive  results. 

A.  UNIFORM  DISTRIBUTION  APPROXIMATION 

By  altering  the  search  scenario  modeled  in  SUBTRAN,  a  uniform 
distribution  of  the  times  to  detection  should  be  achieved.  This  contention  is 
obvious  from  the  fact  that  at  the  start  of  each  simulation  replication  the  target 
is  uniformly  positioned  in  the  search  area.  With  the  target’s  speed  set  at  zero, 
the  time  it  takes  the  searcher  to  detect  the  target  on  one  complete  search  of  the 
area  should  be  uniformly  distributed  on  the  open  interval  (0,T).  Where  T  is  the 
time  it  takes  the  searcher  to  make  one  complete  search  of  the  area. 

This  uniform  distribution  approximation  was  tested  on  SUBTRAN  using  four 
different  sets  of  input.  The  input  values  are  summarized  in  Table  IV.  1.  Figures 
4.1  through  4.4  graphically  compare  the  estimated  probability  of  non-detection  as 
a  function  of  search  time  with  the  theoretical  probability  of  non-detection  for  the 
four  cases.  The  theoretical  probability  of  non-detection  is  indicated  as  a  dotted 
line.  The  estimated  mean  time  to  detection  (M'lTD),  its  95%  confidence  interval, 
and  the  theoretical  MTTD  are  also  included. 


TABLE  IV.l.  UNIFORM  APPROXIMATION  INPUT  PARAMETERS. 


Case  1 

R 

L 

s. 

s. 

1 

(nm) 

(nm) 

(nm/h) 

(nm/h) 

(h) 

1  1 

10 

100 

0 

10 

48 

2  1 

5 

100 

0 

10 

99 

3  1 

10 

100 

0 

5 

96 

4  1 

10 

225 

0 

10 

264.5 

*  1000  replications  performed 

in  all 

cases. 
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The  results  of  the  uniform  approximation  test  strongly  indicate  proper 
functioning  of  the  simulation  with  regard  to  searcher  movement  and  target 
positioning.  One  note  worth  mentioning  deals  with  the  slight  graphical  deviation 
in  theoretical  and  estimated  probability  of  non-detection  shovm  in  Figure  4.2, 
case  4.  The  graphical  results  indicate  that  detection  in  the  simulated  case  occurs 
sooner  than  detection  in  the  theoretical  case  at  large  values  of  time.  The  reason 
for  this  is  explained  by  the  geometry  of  the  particular  search.  In  the  first  three 
test  cases  the  length  dimension  of  the  search  area  was  evenly  divisible  by  the 
searcher’s  track  width.  As  a  result,  a  complete  search  of  the  area  required  no 
overlap  of  search  efibrt.  However,  in  the  fourth  case  tested,  the  division  of  the 
length  dimension  by  the  track  width  did  not  result  in  a  integer  and  overlap 
occurred  on  the  last  leg.  This  overlap  caused  the  deviation  in  the  CDF.  Target 
detections  that  should  have  theoretically  taken  place  on  the  last  search  leg 
occurred  one  leg  earlier. 

B.  RANDOM  TOUR  MODEL  APPROXIMA-nON 

Shifting  to  a  moving  target  and  a  stationary  searcher,  the  results  of 
SUBTRAN  can  be  compared  directly  with  the  random  tour  model  [Ref.  10:pp.  59- 
62].  The  random  tour  model  was  developed  using  Monte  Carlo  simulation 
techniques  and  applies  to  the  detection  of  a  randomly  moving  target  constrained 
to  a  square  search  area  by  boundary  reflection.  The  model  edso  assumes  that  the 
time  between  target  course  changes  is  exponentially  distributed  with  rate 
parameter  t  and  that  the  searcher  is  fixed  in  the  center  of  the  search  area  with 
a  cookie-cutter  detection  radius.  The  output  of  the  model  estimates  the 
probability  of  non-detection  with  respect  to  time  for  given  values  of  searcher 
detection  range  (R),  target  speed  (V),  search  area  size  (A),  and  the  target  course 
change  rate  parameter  (x).  The  equation  for  the  probability  of  non-detection 
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L  =  100  nm 


Figure  4.1.  Uniform  Approximation  Case  1  >  Probability  of 
Non-Detection. 
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Figure  4.2.  Uniform  Approximation  Case  2  •  Probability  of 
Non-Detection. 
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PROBABIUrr  OF  NON-DETECnON 


Non«Detection. 


R 

L 


=  10  run 
=  225  nm 


Figure  4.4.  Uniform  Approximation  Case  4  •  Probability  of 
Non-Detection. 
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(PND),  given  that  detection  did  not  occur  at  time  zero,  proposed  by  the  random 
toiur  model  is; 

PND(t)  =  exp"^ 

where, 

a  =  [(24.7RV*)  /  t(A  -  x  [1  -  exp(-0.084tA‘’‘A0] 

Four  SUBTRAN  simulation  runs  were  compared  to  the  random  tour  model. 
The  input  parameters  are  summarized  in  Table  IV.2.  Figures  4.5  through  4.8 
graphically  compare  the  estimated  probabihty  of  non-detection  as  a  function  of 
search  time  with  the  random  tour  model  estimate  for  probability  of  non-detection 
for  the  four  cases.  The  random  tour  model  estimate  is  indicated  as  a  dotted  line. 
The  estimated  MTTD,  its  95%  confidence  interval,  and  the  random  tour  model 
(RTM)  estimate  of  MTTD  are  also  included. 


TABLE  IV.2.  RANDOM  TOUR  MODEL  APPROXIMATION  INPUT 
PARAMETERS. 


Case  1 

1 

R 

(nm) 

L 

(lun) 

s. 

(nm/h) 

T 

(h-') 

(h) 

1  1 

10 

100 

10 

1 

1000 

2  1 

5 

100 

10 

1 

1000 

3  1 

10 

100 

5 

1 

1000 

4  1 

10 

100 

10 

0.6 

1000 

*  1000  replications  were 

performed  in  all  cases. 

This  test  shows  the  close  comparison  between  the  results  produced  by 
SUBTRAN  and  those  predicted  by  the  random  tour  model.  In  all  cases  tested, 
SUBTRAN  produced  a  slightly  shorter  time  to  detection  than  the  random  tour 
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C.5.  Random  Tour  Model  Approximation  Case  1  •  Probability  of 
Non>Detection. 
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Figure  4.6.  Random  Tour  Model  Approximation  Case  2  •  Probability  of 
Non>Detection. 
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PROBABILITY  OF  NON-DETECTION  3  PROBABIUTY  OF  NON-DETECTION 
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.7.  Random  Tour  Model  Approximation  Case  3  •  Probability  of 
Non-Detection. 
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Figure  4.8.  Random  Tour  Model  Approximation  Case  4  •  Probability  of 
Non-Detection. 
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model.  These  results  strongly  validate  the  modeling  of  target  movement,  target 
boundary  reflection,  and  computation  of  the  occurrence  of  steady-state  detection. 


M 


m 
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V.  EMPLOYMENT  OF  SUBTRAN 
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A.  INTRODUCTION 

This  chapter  is  designed  to  demonstrate  the  employment  of  SUBTRAN.  Five 
separate  search  scenarios  involving  plausible  input  parameters  were  simulated 
using  various  combinations  of  target  and  false  transient  occurrence  rates.  The 
output  statistics  of  these  five  test  cases  are  presented  in  graphical  form  to 
emphasize  the  relationship  between  mean  time  to  detection  (MTTD),  mean  time 
between  target  transients  (MTBTT),  and  mean  time  between  false  transients 
(MTBFT).  In  addition,  the  distributions  of  the  times  to  detection  are  analyzed  to 
assess  the  assumption  that  they  belong  to  the  exponential  class  of  distributions. 

Two  additional  search  scenarios,  Excursion  1  and  Excursion  2,  were 
simulated  to  provide  a  further  understanding  of  how  false  transient  occurrences 
affect  MTTD. 

The  numeric  output  statistics  generated  by  SUBTRAN  were  used  to  develop 
the  graphs  displayed  in  this  chapter. 

B.  INPUT  PARAMETERS 

The  input  parameters  that  remained  common  among  all  five  test  case  and 
the  two  excursions  are: 

-  Steady-state  sotmd  propagation  loss  curve  displayed  in  Figure  5.1. 

-  Transient  sound  propagation  loss  curve  equal  to  the  steady-state  sound 

propagation  loss  curve  displayed  in  Figure  5.1. 

-  Transient  sound  figure  of  merit  set  at  90  dB. 

-  Target’s  course  cumulative  distribution  curve  displayed  in  Figure  5.2. 
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•  Target’s  time  to  course  change  cumulative  distribution  curve  displayed  in 

Figure  5.3. 

-  Lambda-Sigma  Jiimp  model  parameters.  Lambda  set  at  3  hour'*;  Sigma  at  6 

dB. 

-  Track  width  flag  set  at  0. 

-  Number  of  simulation  replications  set  at  1000. 

The  input  parameters  that  varied  over  the  five  test  cases  and  two 
excursions  were;  the  mean  steady-state  detection  range  (R),  the  length  dimension 
of  the  search  area  (L),  the  target’s  speed  (S,),  the  searcher’s  speed  (S,),  and  the 
maximum  search  time  (T^^).  The  values  of  these  parameters  are  summarized  in 
Table  V.l. 

The  target  and  false  transient  occurrence  rates,  converted  to  MTBTT  and 
MTBFT  respectively,  used  in  each  test  case  and  excursion  are: 

-  MTBTT’s  of  1000,  200,  100,  75,  50,  10,  6  (in  hours). 

-  M'TBFT’s  of  400,  75,  10,  5,  1  (in  hours). 

In  addition,  simulation  runs  were  made  in  which  no  target  and/or  false 
transient  sounds  were  allowed. 

C.  OUTPUT  STATIS-nCS  AND  ANALYSIS 

The  MTTDs  and  associated  95%  confidence  intervals  for  the  five  test  cases 
and  two  excursions  are  tabulated  in  Tables  V.2  through  V.8. 

Three  dimensional  graphical  displays  of  the  relationship  between  MTTD, 
MTBIT,  and  MTBFT  for  the  five  test  cases  and  two  excursions  are  contained  in 
Figures  5.4  through  5.8,  5.14,  5.16. 

Graphical  analyses  of  the  probability  of  non-detection  (PND)  for  the  five  test 
cases  and  two  excursions  are  displayed  in  Figures  5.9  through  5.13,  5.15,  5.17. 


Figure  6.1.  Steady-State  and  Transient  Sound  Propagation  Loss  Curve. 
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CUMMUUTIVE  PROBABIUTY 


COURSE  (OEG) 


Figure  6.2.  Target’s  Course  Cumulative  Distribution  Curve. 


Figure  6.3.  Target’s  Time  to  Course  Change  Cumulative  Distribution 
Curve. 


TABLE  y.l.  EMPLOYMENT  INPUT  PARAMETERS. 


Test  Case  1 

R 

L 

s» 

s. 

1 

(nm) 

(nm) 

(nm/h) 

(nm/h) 

(h) 

1  1 

100 

1000 

2  1 

2.5 

100 

1000 

3  1 

2.5 

100 

1000 

4  1 

2.5 

100 

2000 

5  1 

2.5 

150 

4000 

EX  1  1 

10.0 

100 

0.25 

10 

1000 

EX  2  1 

5.0 

100 

5.0 

10 

1000 

*  EIX  2  involves  deceptive  false  transient  sounds. 

1.  Test  Cases  1  Through  6 

Test  cases  1  through  5  involve  variations  in  R,  L,  St,  S„and  T„„.  The 
variations  of  R,  L,  S.,  and  S,  were  chosen  to  encompass  plausible  search 
scenarios.  T„„  was  chosen,  based  on  the  search  scenario  involved,  to  minimize 
the  number  of  simulation  runs  in  which  steady-state  detection  did  not  occur, 
a.  Relationship  Between  MTTD,  MTBTT,  and  MTBFT 
Figfures  5.4  through  5.8  indicate  the  following  three  points; 

-  MTTD  decreases  as  MTBTT  decreases  throughout  the  range  of  MTBTT’s 
tested. 

-  The  "decrease"  trend  is  present  in  MTTD,  as  a  function  of  MTBTT,  for  all 
values  of  MTBFT  tested. 

-  With  respect  to  MTBFT,  MTTD  is  relatively  constant  for  a  given  MTBTT 
unless  a  high  false  transient  occurrence  rate  is  present  (i.  e.  MTBFT  <  5 
hours). 


These  three  points  hif^ilight  the  general  results  indicated  by  Test 
Cases  1  through  5.  SpeciiicaUy,  target  transient  sounds  occurring  at  relatively 
moderate  rates  aid  the  searching  submarine  in  the  steady-state  search  effort. 
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while  false  transient  sounds  only  hamper  the  search  effort  when  their  occurrence 
rate  is  high. 

Figures  5.9  through  5.13  further  emphasize  the  effects  of  target  and 
false  transient  sounds  on  the  searcher’s  steady-state  search  effort.  Three  lines 
are  plotted  on  each  graph  representing  the  PND  for  the  three  scenarios:  no 
target  or  false  transient  sounds  present,  only  target  transient  sounds  present 
with  a  MTBTT  of  5  hours,  and  both  target  and  false  transient  sounds  present 
with  a  MTBTT  of  5  hours  and  a  MTBFT  of  1  hour.  The  improvement  in  the 
time  to  steady-state  detection  is  indicated  by  comparing  the  line  associated  with 
no  transient  sounds  and  the  line  associated  with  only  target  transient  sounds. 
The  degradation  in  the  time  to  steady-state  detection  caused  by  the  introduction 
of  false  transients  is  indicated  by  comparing  the  line  associated  with  only  target 
transient  sounds  and  the  line  with  both  target  and  false  transient  sounds.  In  all 
of  the  test  cases  the  magnitude  of  the  decrease  in  the  time  to  steady-state 
detection  realized  by  the  introduction  of  target  transient  sounds  with  a  MTBTT 
of  5  hours  is  greater  than  the  magnitude  of  the  increase  in  the  time  to  steady- 
state  detection  realized  by  adding  false  transient  sounds  with  a  MTBFT  of  1 
hour. 

b.  Distribution  of  the  Times  to  Detection 

Figures  5.9  through  5.13  can  also  be  used  to  graphically  assess  the 
exponential  "fit"  of  the  estimated  PNDs.  Theoretical  PND  curves,  assuming  the 
times  to  detection  are  exponentially  distributed  with  a  mean  equal  to  the 
estimated  mean  time  to  detection,  are  plotted  as  a  dashed  lines  for  comparison. 
It  is  apparent  that  for  the  five  test  cases,  PND  can  be  closely  approximated  by 
an  exponential  distribution  with  a  mean  equal  to  the  mean  time  to  detection. 


Additional  support,  in  favor  of  the  exponential  approximation  of  the 
PND,  was  provided  by  examining  the  output  statistics  generated  by  SUBTRAN. 
The  estimated  mean  of  the  times  to  detection  and  its  associated  95%  confidence 
interval,  the  estimated  standard  deviation  of  the  times  to  detection,  and  the 
mean  time  to  detection  calculated  using  linear  regression  of  the  times  to 
detection  (see  Chapter  III.C.2  for  explanation)  were  numerically  compared  for 
agreement.  If  the  underlying  distribution  is  truly  exponential,  all  three  vadues 
should  be  approximately  equal.  Of  the  215  simulations  involved  in  the  five  test 
cases,  205  simulations  resulted  in  the  difference  between  the  estimated  mean  of 
the  times  to  detection  and  the  estimated  standard  deviation  to  be  within  the 
95%  confidence  interval  associated  with  the  mean,  while  206  simulations 
indicated  similar  agreement  between  the  estimated  mean  of  the  times  to 
detection  and  the  mean  time  to  detection  calculated  using  linear  regression.  The 
R*  confidence  factor  in  all  simulations  was  greater  than  0.98  with  the  meyority  of 
the  values  greater  than  or  equal  to  0.99. 

2.  Excursions  1  and  2 

In  an  effort  to  better  understand  when  false  transients  play  a  major 
role  in  the  steady-state  search  effort,  two  excursion  scenarios  were  simulated. 
The  first  excursion  involved  a  slow,  almost  stationary  (speed  of  0.25  kts/hr), 
target  and  a  searcher  with  a  large  mean  steady-state  detection  range.  The 
second  excursion  involved  "deceptive"  false  transient  sounds. 

In  the  first  excursion,  an  attempt  was  made  to  test  a  search  scenario  in 
which  false  transient  sounds  should  noticeably  degrade  the  search  effort.  As 
presented  in  Chapter  IV,  the  times  to  detection  are  uniformly  distributed  in  a 


MTTD  (HRS) 


TABLE  V.2.  TEST  CASE  1  -  ESTIMATED  MTTD  VALUES  WITH  96%  Cl. 


MTBFT 


1 

5 

10 

75 

400 

NONE 

5 

1 

38±2 

31±1 

29±1 

27±1 

28±1 

29±1 

M 

10 

1 

54±3 

43±2 

41±2 

39±2 

39±2 

42±2 

T 

60 

1 

74±4 

68±4 

66±3 

67±3 

66±3 

63±3 

B 

75 

1 

75±4 

69±4 

71±4 

71±4 

67±4 

70±4 

T 

100 

1 

76±4 

71±4 

74±4 

74±4 

76±4 

75±4 

T 

200 

1 

79±4 

80±4 

75±4 

76±4 

80±4 

75±4 

1000 

1 

84±5 

82±4 

84±5 

80±4 

83±4 

83±4 

NONE 

1 

- 

• 

- 

• 

. 

85±5 

*  All  values  of  MTBFT,  MTBTT,  and  MTTD  are  in  hours 
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MTTD  (HRS) 


Figure  6.5.  Test  Case  2  •  Estimated  MTTD,  M'fBlT,  and  MTBFT 

Relationship. 


TABLE  V.3.  TEST  CASE  2  -  ESTIMATED  MTTD  VALUES  WITH  95%  Cl. 


MTBFT 


1 

5 

10 

75 

400 

NONE 

5 

1 

59±3 

44±2 

42±2 

41±2 

43±2 

42±2 

M 

10 

1 

89±5 

71±3 

69±3 

65±3 

65±3 

69±3 

T 

60 

1 

138±7 

132±7 

124±6 

120±5 

124±6 

122±6 

B 

75 

1 

146±8 

135±6 

135±7 

122±6 

133±7 

128±7 

T 

100 

1 

150±8 

134±7 

138±7 

137±7 

146±7 

135±7 

T 

200 

1 

152±8 

142±7 

142±7 

145±7 

150±8 

146±8 

1000 

1 

156±8 

153±8 

150±8 

158±8 

154±8 

166±8 

NONE 

1 

- 

. 

. 

. 

- 

157a:8 

*  AU  values  of  MTBFT,  MTBTT,  and  MTTD  are  in  hours 


Figure  6.6.  Test  Case  3  •  Estimated  MTTD,  MTBTT,  and  MTBPT 

Relationship. 


TABLE  V.4.  TEST  CASE  3  -  ESTIMATED  MTTD  VALUES  WITH  96%  Cl. 


MTBFT 


1 

5 

10 

75 

400 

NONE 

5 

i 

37±2 

28±1 

28±1 

29±1 

28±1 

27+1 

M 

10 

1 

57±3 

48±3 

47±2 

44+2 

44±2 

43±2 

T 

50 

1 

102±5 

94±5 

98±5 

95±5 

90±5 

90±4 

B 

76 

1 

105±6 

108±6 

102±5 

1081:6 

101±5 

103±6 

T 

100 

1 

106±6 

108±6 

108+5 

112±6 

107±5 

103±5 

T 

200 

1 

119±6 

120±6 

118±6 

122±6 

116±6 

114±6 

1000 

1 

124±7 

127±7 

126±6 

126±7 

121±6 

130±7 

NONE 

1 

. 

- 

- 

127±7 

*  All  values  of  MTBFT,  MTBTT,  and  MTTD  are  in  hours 
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MTTD  (HRS) 


Relationsliip. 


TABLE  V.6.  TEST  CASE  4  •  ESTIMATED  MTTD  VALUES  WITH  95%  Cl. 


MTBFT 


1 

5 

10 

75 

400 

NONE 

5 

1 

95±5 

53±2 

51±2 

47±2 

50±2 

49±2 

M 

10 

1 

150±7 

86±4 

84±4 

81±3 

79±3 

81±3 

T 

60 

1 

243±13 

198±10 

186±9 

185±9 

174±9 

192±10 

B 

75 

1 

264i;13 

224±11 

206±10 

212±11 

201±10 

221ill 

T 

100 

1 

268±14 

234±13 

211±10 

219±11 

218±12 

230±11 

T 

200 

1 

281±15 

243±12 

232±12 

252±13 

243±13 

238±13 

1000 

1 

289±15 

276±13 

2601 14 

269±15 

288±15 

271±16 

NONE 

1 

- 

• 

- 

- 

- 

273±14 

*  All  values  of  MTBFT,  MTBTT,  and  MTTD  are  in  hours 
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MTTD  (HRS) 


2.5  nm 
150  nm 
3  nm/hr 
5  nm/hr 
4000  hr 


Figure  5.8.  Test  Case  5  -  Estimated  MTTD,  MTBTT,  and  MTBFT 

Relationship. 


TABLE  V.6.  TEST  CASE  5  -  ESTIMATED  MTTD  VALUES  WITH  96%  CL 


MTBFT 


1 

5 

10 

75 

400 

NONE 

5 

1 

152±7 

108±4 

108a:4 

111±5 

111±5 

113±5 

M 

10 

1 

266±14 

169±8 

163±8 

165±7 

169±7 

173±7 

T 

50 

1 

494±26 

416±21 

416±20 

403±19 

389±18 

392±18 

B 

75 

1 

548±29 

470±23 

446±21 

441±22 

432±22 

477±23 

T 

100 

1 

553±29 

472±24 

473±24 

474±24 

470±25 

481±24 

T 

200 

1 

597±32 

552±27 

542±26 

583±30 

579±29 

657±28 

1000 

I 

688±35 

582±29 

579±30 

637±32 

614±31 

601±32 

NONE 

1 

• 

- 

* 

. 

604±31 

*  All  values  of  MTBFT,  MTBTT,  and  MTTD  are  in  hours 
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PR0BAB1U7Y  OF  NON-DETECTION  PROBABILITY  OF  NON-DETECTION 


1  •  no  transients 

2  '  target  transients  only 


Figure  5.9.  Test  Case  1  •  Probability  of  Non-Detection. 


1  -  no  transients 


2  •  target  transients  only 
(MTBTT  =  5  hrs) 


Figure  6.10.  Test  Case  2  •  Probability  of  Non-Detection. 


1  -  no  transients 

2  -  target  transients  only 


Figure  5.H.  Test  Case  3  -  Probability  of  Non-Detection. 


1  -  no  transients 

2  -  target  transients  only 


Figure  6.12.  Test  Case  4  -  Probability  of  Non-Detection. 
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1  •  no  transients 

2  -  target  transients  only 


Figure  5.13.  Test  Case  5  •  Probability  of  Non>Detection. 


scenario  involving  a  stationary  target,  a  systematically  moving  searcher,  and  no 
transient  sounds.  However,  by  introducing  target  and  false  transient  sounds  to 
the  scenario,  the  searcher  is  forced  away  from  his  systematic  search  by  transient 
occurrence  investigations.  In  essence,  the  search  effort  is  degraded  to  a  random 
search  of  the  area. 

Comparing  the  mean  time  to  detection  in  the  uniform  distribution  case 
to  the  mean  time  to  detection  in  the  random  search  case  we  have; 

MTTD„,^orm=  A/4RVe  MTTD«,^«,m=  A  /  2RV, 

where, 

A  =  area  of  search 


R  =  radius  of  detection 
Vg  =  effective  search  speed 
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Since  the  mean  time  to  detection  in  the  uniform  distribution  case  is  one 
half  that  in  the  random  search  case,  the  degradation  in  MTTD  experienced  with 
the  introduction  of  false  transient  sounds  should  be  greater  than  the  degradation 
that  occurred  in  the  five  test  cases.  Figure  5.14  indicates  this  contention.  In 
addition,  Figure  5.15,  indicates  the  transition  from  a  xmiform  search  to  a  random 
search. 

In  Figure  5.15,  the  curve  associated  with  no  false  transient  sounds  has 
deviated  from  the  exponential  fit  as  is  indicated  by  the  bowing  of  the  curve  on  a 
natural  log  scale.  Finally,  Figure  5.15  also  refutes  the  trend,  established  in  the 
five  test  cases,  that  the  gains  achieved  by  target  transient  sounds  were  greater 
than  the  losses  experienced  from  the  introduction  of  false  transient  sounds  at 
nearly  equal  rates  of  occurrence.  The  curve  associated  with  both  target  and  false 
transient  sounds  shifts  back  to  the  position  of  the  curve  associated  with  no 
transients  and  at  higher  time  values  surpasses  it. 

Excursion  2  deals  with  "deceptive"  false  transients.  The  main  program  of 
SUBTRAN  was  altered  slightly  to  achieve  a  method  for  drawing  the  searcher 
away  from  the  target  in  the  event  of  a  false  transient  sound  emission.  In 
Excursion  2,  at  the  time  of  a  false  transient  occurrence,  the  position  of  false 
transient  occurrence  (X^Yf)  was  set  at: 

Xf  =  L  -  X. 

Yr  =  L  -  Y, 

Using  this  method,  the  searcher  was  directed  away  from  the  target  in 
most  cases.  The  results  of  this  excursion  are  graphically  displayed  in  Figures 
5.16  and  5.17.  The  figures  clearly  indicate  the  increased  effect  that  false 
transient  sounds  have  on  the  search  effort.  In  particular,  it  is  interesting  to  note 
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in  Figure  5.17  that  the  ctirve  associated  with  both  target  and  false  transient 
sounds  has  completely  shifted  past  the  curve  associated  with  no  transient 
sounds. 
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MTTD  (HRS) 


TABLE  V.7.  EXCURSION  1  -  ESTIMATED  MTTD  VALUES  WITH  96%  Cl. 


MTBFT 


1 

5 

10 

75 

400 

NONE 

5 

1 

21±1 

15±1 

14±1 

13±1 

13±1 

13±1 

M 

10 

1 

25±2 

19±1 

18±1 

16±1 

16±1 

15±1 

T 

50 

( 

28±2 

25±1 

22±1 

19±1 

18±1 

18±1 

B 

75 

1 

28±2 

26±1 

24±1 

19±1 

19±1 

19±1 

T 

100 

1 

33±2 

27±2 

24±1 

19±1 

20±1 

19±2 

T 

200 

1 

31±2 

26±1 

24i:l 

19±1 

19±1 

19±1 

1000 

1 

31±2 

27±2 

241:1 

20±1 

20±1 

19±1 

NONE 

1 

- 

. 

. 

19±1 

*  All  values  of  MTBFT,  MTBTT,  and  MTTD  are  in  hours 
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PROBABIUTY  OF  NON-DETECTION 


1  '  no  transients 

2  '  target  transients  only 

(MTBTT  =  5  hrs) 


Figure  6.15.  Excursion  1  -  Probability  of  Non-Detection. 
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MTTD  (HRS) 


(^RS) 


Figure  5.16.  Excursion  2  -  Estimated  MTTD,  MTBTT,  and  MTBFT 

Relationship. 


TABLE  V.8,  EXCURSION  2  -  ESTIMATED  MTTD  VALUES  WITH  96%  Cl. 


MTBFT 


1 

5 

10 

75 

400 

NONE 

5 

1 

125±7 

39+2 

33+1 

27+1 

27±1 

29±1 

M 

10 

1 

173±9 

59+3 

49±3 

41+2 

38±2 

42±2 

T 

50 

1 

221±11 

111±6 

86±4 

68±3 

68±3 

63±3 

B 

75 

1 

224±12 

112±6 

89+5 

70±3 

71+4 

70±4 

T 

100 

1 

233±12 

120±6 

90+5 

75±4 

72±4 

75±4 

T 

200 

1 

226+11 

119+6 

96+5 

83+4 

75±4 

75±4 

1000 

1 

230±12 

127±7 

109+6 

82±4 

78±4 

83±4 

NONE 

1 

- 

. 

- 

- 

- 

85±6 

*  All  values  of  MTBFT,  MTBTT,  and  MTTD  are  in  hours 
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PROBABIUrr  OF  NON-DETECnON 


1  -  no  transients 

2  •  target  transients  only 


Figure  5.17.  Excursion  2  •  Probability  of  Non-Detection. 
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VL  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  WORK 


Although  this  thesis  does  not  provide  an  all-inclusive  answer  to  the  question 
of  whether  or  not  transient  detection  opportuxiities  can  be  utilized  by  a  seeu'ching 
submarine  to  increase  its  chance  of  gaining  steady-state  detection  on  a  target 
submarine,  it  does  establish  a  framework,  through  the  simulation  SUBTRAN, 
that  may  be  useful  in  investigating  how  transient  sounds  affect  the  steady-state 
search  effort  for  a  given  search  scenario. 

The  modeling  in  SUBTRAN  utilizes  established  methods.  The  target 
movement  and  transient  occurrence  models  allow  the  user  flexibility  in  defining 
the  target’s  characteristics.  The  searcher  movement  model  is  restrictive  in  the 
sense  that  the  searcher  has  no  movement  options  in  the  search  pattern  or 
transient  prosecution.  The  false  transient  occurrence  model  is  designed  to 
supplement  the  search  scenario  with  randomly  occurring  ambient  sounds  that 
could  be  classified  as  target  transient  sounds  by  the  searcher. 

Major  features  of  SUBTRAN  that  support  it  as  a  viable  method  for 
simulation  of  the  transient  environment  search  scenario  presented  in  this  thesis 
are; 

-  Considerable  flexibility  is  available  to  the  user  in  selection  of  search  scenario 
input  parameters. 

-  The  simulation  is  designed  for  microcomputer  implementation  removing  the 
"mainframe"  requirement  common  among  most  submarine  versus  submarine 
simulation  programs. 

-  The  simulation  search  scenario  is  easy  to  set-up  using  the  input  data  file 
provided. 
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-  The  output  provided  by  SUBTRAN  is  concise  and  clearly  highlights 
important  results. 

Improvements  or  additions  SUBTRAN  might  include: 

-  Modeling  target  counter-detection  of  the  searcher. 

-  Incorporating  different  types  of  target  transients  to  include  various  duration 
length  and  FOMs. 

-  Investigating  different  searcher  movement  patterns. 

-  Including  FOM  target  aspect  and  speed  dependencies. 

-  Modeling  target  and  searcher  movement  to  include  speed  changes. 
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APPENDIX  A 


SUBTRAN  USERS  GUIDE 

A.  SYSTEM  REQUIREMENTS 

SUBTRAN  requires  both  computer  memory  and  disk  space.  The  computer 
memory  required  is  126  Kbjrtes.  The  disk  space,  necessary  for  the  two  output 
files,  varies  with  respect  to  the  number  of  simulation  replications  performed.  The 
disk  space  needed  for  a  simulation  run  of  6000  is  approximately  101  Kbytes, 
while  the  disk  space  needed  for  a  simulation  run  of  1000  is  approximately  27 
Kbytes. 

B.  INPUT  DATA  FILE 

SUBTRAN  requires  input  parameters  to  be  entered  via  an  input  data  file. 
The  input  data  file  is  specially  formatted  to  aid  the  user  in  proper  positioning  of 
the  input  parameters.  A  sample  of  the  input  data  file  is  contained  in  Appendix 

C. 

SUBTRAN  was  not  coded  to  check  for  improper  input  parameters.  Careful 
attention  is  required  when  filling  out  the  input  data  file  to  ensure  input 
parameters  are  within  specified  range  limits. 

The  input  data  file  is  broken  up  into  15  data  fields.  A  description  of  the 
data  requirements  and  allowed  ranges  follows; 

1.  Data  Field  1 

The  input  parameters  contained  in  data  field  1  are:  random  number 
seed  (DSEED),  length  dimension  of  the  search  area  (L),  number  of  simulation 
replications  (N),  and  the  maximum  search  time  (T„„). 


DSEED  -  double  precision  variable  range,  [1,2147483646] 

L  -  real  variable  /  range,  (0,999] 

N  -  real  variable  /  range,  [1,5000] 

-  real  variable  /  range,  [0,9999] 

2.  Data  Field  2 

The  input  parameter  contained  in  data  Held  2  is  the  number  of  points 
to  be  entered  for  the  steady-state  propagation  loss  curve  (NPp,,). 

NPp.,  -  integer  variable  /  range,  [0,50] 

3.  Data  Field  3 

The  input  parameters  contained  in  data  field  3  are  the  steady-state 
propagation  loss  curve  range  values  in  increasing  order  (RN^)  and  the 
corresponding  propagation  loss  values  (Pp.,).  Due  to  the  code  structure  for 
calculating  the  propagation  loss  curve,  the  values  for  any  two  sequential  range 
inputs  must  be  different.  Also,  the  last  propagation  loss  value  must  be  the 
highest  of  all  propagation  loss  values. 

RNp^  -  real  array  of  size  NP^^  /  range,  [0,999] 

Pp.  -  real  array  of  size  NPp.  /  range,  [0,999] 

4.  Data  Fields  4  and  5 

The  input  parameters  contained  in  data  fields  4  and  5  are  of  the  same 
structure  as  the  parameters  contained  in  data  fields  2  and  3  respectively. 
Instead  of  steady-state  propagation  loss  curve,  the  input  corresponds  to  the 
transient  propagation  loss  curve.  The  variables  in  data  fields  4  and  5 
corresponding  to  those  in  data  fields  2  and  3  are; 
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NP,., 

replaces 

NPp., 

RNp. 

replaces 

RN,. 

Pp. 

replaces 

Pp- 

5.  Data  Field  6 

The  input  parameters  contained  in  data  field  6  are  the  searcher’s  speed 
(S.)  and  the  target’s  speed  (T.).  The  speed  parameters  must  never  be  set  to  zero. 

S.  -  real  variable  /  range,  (0,99] 

S,  -  real  variable  /  range,  (0,99] 

6.  Data  Field  7 

The  input  parameter  contained  in  data  field  7  is  the  number  of  points 
to  be  entered  for  the  target’s  course  distribution  (NP,.). 

NP,  -  integer  variable  /  range,  [0,50] 

7.  Data  Field  8 

The  input  parameters  contained  in  data  field  8  are  the  target’s  course 
distribution  values  (C^)  in  increasing  order  and  the  corresponding  cumulative 
probability  (P,,).  The  first  value  of  the  cumulative  probability  must  be  0  and  the 
last  value  of  cumulative  probability  must  be  1. 

C,,  -  real  array  of  size  NP,.  /  range,  [0,360] 

Pj  -  real  array  of  size  P,,  /  range,  [0,1] 

8.  Data  Fields  9  and  10 

The  input  parameters  contained  in  data  fields  9  and  10  are  of  the  same 
structure  as  those  contained  in  data  fields  7  and  8  respectively.  Instead  of 
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target’s  course  distribution,  the  input  corresponds  to  the  distribution  of  the  time 
between  target  course  changes.  The  variables  in  data  fields  9  and  10  that 
correspond  to  those  in  data  field  7  and  8  are: 

NPfe  replaces  NP, 

replaces  C„  range  for  T,.  is  10,99] 

P^  replaces  P,, 

9.  Data  Field  11 

The  input  parameters  contained  in  data  field  11  are  the  transient 
occurrence  flag  (F^)  and  the  transient  occurrence  rate  parameter  (RP„).  Transient 
occurrence  is  allowed  when  the  transient  occurrence  flag  is  set  to  a  value  of  1. 
Transients  are  not  allowed  when  the  flag  is  set  to  a  value  of  0.  The  value  of  the 
transient  occurrence  rate  parameter  must  never  be  0. 

F„  -  integer  variable  /  range,  0  or  1 

RP^  -  real  variable  /  range,  [.000001,99] 

10.  Data  Field  12 

The  input  parameters  contained  in  data  field  12  are  of  the  same 
structure  as  those  input  parameters  contained  in  data  field  11.  Instead  of 
transient  occurrence,  the  input  corresponds  to  false  transient  occurrence.  The 
variables  in  data  field  12  that  correspond  to  those  in  input  data  field  11  are; 

Fa,  replaces  F^ 

RPft,  replaces  RP^ 
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11.  Data  Field  13 


The  input  parameters  contained  in  input  data  field  13  are  the  Lambda- 
Sigma  Jump  error  process  rate  parameter  (RP,^)  and  the  Lambda-Sigma  Jump 
error  process  standard  deviation  (SD^,).  The  rate  pvameter  must  never  be  0. 

RP,^  -  real  variable  /  range,  t-000001,99] 

SD,,j  -  real  variable  /  range,  [0,99] 

12.  Data  Field  14 

The  input  parameters  contained  in  input  data  field  14  are  the  steady- 
state  figure-of-merit  (FOM„)  and  the  transient  figure-of-merit  (FOM^,). 

FOM,.  -  real  variable  /  range,  [0,999] 

FOMt,  -  real  variable  /  range,  [0,999] 

13.  Data  Field  15 

The  input  parameters  contained  in  input  data  field  15  are  the  track 
width  flag  (F,,)  and  the  track  width  (TW).  The  input  track  width  is  used  in  the 
simulation  if  the  value  of  the  track  width  flag  is  set  at  1.  The  track  width  is  set 
to  the  mean  steady-state  detection  range  if  the  value  of  the  track  width  flag  is 
set  at  0. 

Ft„  -  integer  variable  /  range,  0  or  1 
TW  -  real  variable  /  range,  [0,999] 

C.  PROGRAM  EXECUTION 

SUBTRAN  is  executed  by  running  the  SUBTRAN.EXE  file.  The  user  is 
prompted  to  enter  the  name  of  the  input  data  file.  The  name  of  the  file  may  be 
any  proper  filename  and  extension  of  total  length  less  than  or  equal  to  12 
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characters.  After  the  input  data  file  name  is  entered,  the  input  parameters  are 
printed  at  the  user’s  terminal  for  verification.  Simulation  execution  commences 
when  the  input  parameters  have  been  verified.  Status  of  the  simulation  is 
printed  at  the  user’s  terminal  every  time  a  tenth  of  the  simulation  is  completed. 
When  all  replications  have  been  run,  output  statistics  are  printed  at  the  user’s 
terminal  and  to  the  output  data  files. 

D.  OUTPUT  DATA  FILES 

Two  output  data  files  are  generated  by  SUBTRAN.  The  first  data  file, 
SUBTRAN.OUT,  contains  the  input  parameters  and  output  statistics.  The  second 
data  file,  TIME.OUT,  contains  the  times  to  detection  ordered  from  shortest  to 
longest.  Times  to  detection  with  a  value  of  zero  are  not  included.  In  addition, 
simulation  replications  that  did  not  result  in  steady-state  detection  before  T„„, 
have  time  to  detection  set  at  T„„. 
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APPENDIX  B 


MAJOR  VARIABLES  USED  IN  SUBTRAN 

The  major  variables  used  in  SUBTRAN  are  provided  in  an  alphabetical 
hsting. 

Cj  =  Course  array  of  the  targets  course  cumulative  distribution  curve 
in  degrees. 

C,  =  Course  the  target  will  be  on  after  the  next  boundary  reflection 
in  degrees. 

C.  =  Course  of  the  searcher  in  degrees. 

Ct  =  Course  of  the  target  in  degrees. 

DSEED  =  Random  number  seed. 

E, ^  =  Signal  excess  fluctuation  level  in  dB. 

Fft,  =  False  transient  occurrence  flag. 

Fer  =  Transient  occurrence  flag. 

F, ^  =  Track  width  calculation  flag. 

FOM„  =  Steady-state  figure-of-merit  in  dB. 

FOMt,  =  Transient  figure-of-merit  in  dB. 

L  =  Length  dimension  of  the  search  area  in  nautical  miles. 

N  =  Number  of  simulation  repUcations. 

Nj  =  Number  of  repUcations  in  which  steady-state  detections  occurred. 

N„  =  Number  of  replications  in  which  no  steady-state  detection 

occurred  before  T„„. 

N,  =  Number  of  repUcations  in  which  steady-state  detection  occurred 
at  time  zero. 

NPj  =  Number  of  points  entered  for  the  target’s  course  cumulative 
distribution  curve. 


NPjm  =  Number  of  points  entered  for  the  steady-state  propagation  loss 
curve. 

NPpt,  =  Number  of  points  entered  for  the  transient  propagation  loss 
curve. 

NP^  s  Number  of  points  entered  for  the  target’s  time  to  course  change 

cumulative  distribution  curve. 

Pe  =:  Probability  array  for  the  target’s  course  oimiilative  distribution 

curve. 

Pp.  =  Propagation  loss  array  for  the  steady-state  propagation  loss 

curve  in  dB. 

Ppt,  =  Propagation  loss  array  for  the  transient  propagation  loss  curve 

in  dB. 

Pt,  =  Probability  array  for  the  target’s  time  to  course  change 

cumulative  distribution  curve. 

Rj  =  Range  array  for  transient  detection  ranges  in  nautical  miles. 

R, ^  =  Steady-state  detection  range  in  nautical  miles. 

RNp^  =  Range  array  for  the  steady-state  propagation  loss  curve  in 

nautical  miles. 

RNpt,  =  Range  array  for  the  transient  propagation  loss  curve  in  nautical 

miles. 

RPft,  =  Rate  parameter  for  the  false  transient  occurrence  distribution  in 

hours’. 

RPi^  =  Rate  parameter  for  the  Lambda-Sigma  Jump  process  in  hours  ’. 

RP^  =  Rate  parameter  for  the  transient  occurrence  distribution  in 

hours’. 

S,  =  Speed  of  the  searcher  in  nautical  miles  per  hour. 

St  =  Speed  of  the  target  in  nautical  miles  per  hour. 

SD,,j  =  Standard  deviation  for  the  Lambda-Sigma  Jump  process  in  dB. 

T, ^  =  Time  to  the  next  signal  excess  fluctuation  in  hours. 

Tp,„  =  The  maximum  search  time  for  each  replication  in  hours. 

T,  =  Time  until  the  next  target  boundary  reflection  in  hours. 
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Tte  =  Time  array  for  the  target’s  time  to  coixrse  change  cumulative 
distribution  curve  in  hours. 

TC.  =  Time  to  the  next  searcher  course  change  in  hours. 

TC,  =  Time  to  the  next  target  course  change  in  hours. 

TTf  =  Time  to  the  next  false  transient  occurrence  in  hours. 

TTt  =  Time  to  the  next  transient  occurrence  in  hours. 

TW  =  Searcher’s  track  width  in  nautical  miles. 

Xf  =  X  component  of  the  false  transient  occurrence  position. 

X,  =  X  component  of  the  searcher’s  position. 

Xt  =  X  component  of  the  target’s  position. 

Yf  =  Y  component  of  the  false  transient  occurrence  position. 

Y,  =  Y  component  of  the  searcher’s  position. 

Y,  =  Y  component  of  the  target’s  position. 
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W  *  *  CO  * 


APPENDIX  C 


EXAMPLE  SUBTRAN  INPUT  DATA  FILE 


cccccccccccccccccccccccccc 
c 

C  The  following  is  the  input  file  for  SUBTRAN. 

C  All  values  should  be  right  Justified  within 

C  their  designated  field.  Any  line  with  a  ’** 

C  or  a  'C  in  it  must  not  be  deleted,  success- 

C  ful  reading  of  the  input  data  depends  on 

C  it.  If  continued  errors  occur  when  the  input 

C  is  read  by  the  main  program,  consult  the 

C  program  documentation  for  the  proper  setup 

C  of  this  file. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCC 
* 

1. 


* 

Columns 

Input 

♦ 

* 

t 

3-15 

DSEED 

[  d] 

* 

18-30 

L 

[  1] 

* 

* 

33-45 

N 

t  n] 

* 

48-60 

Tmax 

[  t] 

t 

*  * 

* 

*  *  *  * 

**ddddddddddddd**l  1 111111111 1  l**nnnnnnnnnnnnn**ttttttttttttt1^* 

2855 

100 

1000 

1000 

2.  t  t 

♦ 

t  t  t 

*  * 

* 

Columns 

I  nput 

* 

* 

* 

* 

3-15 

NPpss 

[  n3 

t 

t  * 

* 

♦  ♦  *  * 

^  * 

* 

**nnnnnnnnnnnnn**  [ NPpss  cannot  exceed  50] 
4 

* 

3.  *  * 

* 

*  *  *  * 

*  * 

* 

* 

Columns 

I  nput 

t 

* 

* 

* 

3-15 

RNpss  CNPp) 

[  r] 

* 

* 

18-30 

Ppss (NPp) 

[p] 

* 

*  * 

* 

t  t  t  t 

t  ♦ 

♦ 

♦♦rrrrrrrrrrrrr^^ppppppppppppp** 
0  46 


oononooooooo 


* 

4. 


*  * 


♦  * 


*  *  *  * 
Columns 


*  *  *  * 


»  *  *  * 
Input 


*  *  * 


♦♦nnnnnnnnnnnnn** 

2 

t 

5.  t  t  * 


3-15  NPptr  In] 

[ NPptr  cannot  exceed  50] 


♦ 

* 

* 

* 

* 

* 


*  * 


*  *  ♦  * 
Columns 

3-15 
18-30 
*  *  *  * 


*  t  t 


Input 


*  »  * 


RNptr(NPp)  [  r] 
PptrtNPp)  [p] 
in  %  t  * 


**rrrrrrrrrrrrr**ppppppppppppp** 


* 

6. 


* 

♦ 

* 

* 


*  * 


*  t 


0 

50 

*  *  •*  * 
Columns 

3-15 
18-30 
t  t  t  t 


50 

120 


t  t  t  t 


*  *  * 
Input 


St 
Ss 
*  * 


**tttttttttttttt1;ssssssssssssstt 


10 


7. 


t 

* 

* 

* 

* 


*  * 


t  * 


t  t  t  t 
Columns 

3-15 
*  *  *  t 


t  if  t  t 


t  t  t  t 


t  t  t 
I  nput 

NPc 
in  if  in 


%  % 


%  % 


in  %  t 


♦*nnnnnnnnnnnnn*t 

6 


C  NPc  cannot  exceed  50] 


* 

* 

* 

* 

* 

% 


in  %  in  in 


[  t] 
[  s] 
t  * 


t  t  %  in 


[  n] 
*  % 


8 .  *  * 

if.ini.inininin% 

t  *  *  t  * 

%  in  in 

i 

* 

Columns 

I  nput 

% 

* 

in 

3-15 

Cc(NPc) 

[  c] 

in 

% 

18-30 

Pc (NPc) 

[p] 

in 

in  in 

%%%%%%%% 

in  in  %  in  % 

in  in  in 

* 

♦*ccccccccccccc**pppppppppppppt* 

0 

0 

20 

.  2 

90 

.  5 

180 

.  8 

280 

.  9 

360 

1 

67 


9. 


*  Columns  Input 

♦  -  - 

♦  3-15  NPtc 

**nnnnnnnnnnnnn**  £  NPtc  cannot  exceed 
5 

* 

10.  »**♦****** 


* 

t 

* 

* 

* 

♦ 


*  * 


*  *  *  * 
Columns 

3-15 
18-30 
*  *  *  * 


Input 


*  *  * 


Ttc(NPtc) 
Ptc(NPtc) 
»  *  *  *  * 


**ttttttttttttt**ppppppppppppp;^* 


* 

11. 


* 

* 

% 

% 

% 

% 


1. 

2. 

*  % 


%  % 


0 

42 
89 

43 
04 

»  »  *  * 
Columns 

3-15 
18-30 
%  %  %  % 


%  %  %  % 


t  %  %  % 


0 

2 

4 

6 

1 

*  *  * 
Input 

Ftr 
RPtr 
%  %  % 


%  % 


*  % 


**nnnnnnnnnnnnn**rrrrrrrrrrrrr** 


* 

12, 


* 

* 

* 

* 

* 


*  * 


%  % 


%  t  %  % 
Columns 

3-15 
18-30 
%  t  %  % 


%  %  %  i(. 


%  K,  %  % 


%  %  % 
Input 

Fftr 
RPftr 
%  %  % 


♦♦nnnnnnnnnnnnn**rrrrrrrrrrrrrt* 


13. 


% 

% 

% 

t- 

* 

% 


0 

%%%%%% 

Columns 


*  *  ♦  * 


*  % 


3-15 
18-30 
%  %  %  % 

♦♦rrrrrrrrrrrrr^^sssssssssssss** 
3  6 


*  * 


t  %  % 
Input 


*  t 


*  * 


RPlsJ 
SDlsJ 
%  %  % 


%  t 


*  *  * 
* 
♦ 

£  nl  * 
*  *  * 
50] 


*  *  * 
* 
* 

£  t]  * 
£  p]  * 
*  *  * 


♦  *  » 
* 
* 

£  n]  * 
£  rl  * 
t  t  * 


*  *  * 
t 
* 

£  n]  * 
[  r]  * 
t  t  t 


*  *  * 
* 

E  r]  t 
£  s]  * 
♦  *  * 
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14. 

* 

*  *  ♦  * 

* 

* 

* 

*  *  *  ♦ 

* 

*  t 

* 

Columns 

I  nput 

* 

* 

* 

♦ 

3-15 

FOMss 

1  s] 

* 

♦ 

18-30 

FOMtr 

[  t] 

* 

* 

t 

* 

*  *  * 

* 

* 

* 

t  t  *  * 

* 

* 

*  * 

* 

**sssssssssssss**ttttttttttttt** 

56 

90 

♦ 

15. 

* 

* 

* 

*  *  *  ♦ 

* 

* 

*  *  *  * 

* 

* 

*  * 

* 

* 

Columns 

Input 

* 

* 

* 

3-15 

Ftw 

t  f] 

* 

* 

18-30 

TV 

[  S] 

* 

* 

* 

*  ♦  ♦  * 

* 

* 

* 

*  *  *  * 

* 

* 

*  t 

* 

f  *»sssssssssssss^* 
0  10 
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APPENDIX  D 


EXAMPLE  SUBTRAN  OUTPUT  STATISTICS 


*  SUBTRAN  INPUT  * 


SIMULATION  PARAMETERS  SEARCH  PARAMETERS 


Replications  =  1000 
Random  Seed  =  2855 

TRANSIENT  PARAMETERS 


(1  -  Yes  /  0  -  No) 

Real  Transients: 
Allowed  =  1 

Rate  Par.  =  0.1 

False  Transients: 
Allowed  =  0 

Rate  Par.  =  6.0 


Search  Area  =  10000 
Search  Time  =  1000 
Track  Width  =  calc. 

SIGNAL  ERROR  PARAMETERS 


Rate  Par.  =  3. 00 

Std.  Dev.  =  6.00 

SPEED  PARAMETERS 


S.  Speed  =  10. 00 

T,  Speed  =  5.00 


STEADY-STATE  FIGURE  OF  MERIT  /  PROPAGATION  LOSS  PARAMETERS 


Steady-state  Figure  of  Merit  =  56.00 


Range  Point  : 

0.  00 

Prop. 

Loss  : 

46.  00 

Range  Point  : 

6.  00 

Prop. 

Loss  : 

70.  00 

Range  Point 

14. 00 

Prop . 

Loss  : 

91. 00 

Range  Point  : 

18. 00 

Prop. 

Loss  : 

92.  00 

TRANSIENT  FIGURE  OF  MERIT  /  PROPAGATION  LOSS  PARAMETERS 


Transient  Figure  of  Merit  =  90.00 


Range  Point  ;  0. 00 
Range  Point  ;  50. 00 


Prop.  Loss  :  50.00 

Prop.  Loss  :  120.00 


TARGET  COURSE  PARAMETERS 


Course : 

0. 

Prob: 

0.  0 

T  ime 

to  Change: 

0.  00 

Prob: 

0. 0 

Course : 

20. 

Prob: 

0.2 

Time 

to  Change: 

0. 42 

Prob: 

0.2 

Course: 

90. 

Prob: 

0.5 

Time 

to  Change: 

0.89 

Prob: 

0.4 

Course: 

180. 

Prob: 

0.8 

Time 

to  Change: 

1. 43 

Prob; 

0.  6 

Course : 

280. 

Prob: 

0.  9 

Time 

to  Change: 

2. 04 

Prob: 

1.0 

Course: 

360. 

Prob: 

1.0 

*  * 

*  *  * 

t  t  * 

*  SUBTRAN 

OUTPUT 

■  DATA  * 

*  * 

t  *  * 

*  »  * 

*  *  *  * 

NO.  CF 

REPS. 

WITH  DETECTION  TIME 

:  OF  ZERO  . . 

0 

NO.  OF  REPS.  WITH  DETECTION  TIME  GREATER  THAN  Tnax  ...  0 

CALCULATED  TRACK  WIDTH  VALUE  .  5.00 


TIME  TO  DETECTION  STATISTICS 


Mean 

Std.  Dev. 

1/Lambda 

R-Squared 


644.19  +/-31.13  (95%  Cl) 
598. 44 
602.55 
0.  99 


Prob.  of  E>et.  By  time 


Prob.  of  Det.  By  time 


0.  1 

74.5 

0.2 

154.0 

0. 3 

257.8 

0.  4 

360.  4 

0.  5 

478. 6 

0. 6 

621 . 1 

0. 7 

762. 7 

0. 8 

1054. 4 

0.  9 

1410. 0 

1.  0 

3411. 8 
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APPENDIX  E 


STEADY-STATE  DETECTION  DERIVATION 

As  a  result  of  the  discrete  event  timing  used  in  SUBTRAN,  steady-state 
detection  is  assessed  by  computing  the  time  it  would  take  the  range  between  the 
target  and  the  searcher  to  decrease  to  the  vjilue  of  the  steady-state  detection 
range.  TLe  calculations  used  in  the  steady-state  detection  assessment 
computation  are  outlined  below. 

Defir.e, 

-  Starting  X  position  of  the  target. 

-  Starting  Y  position  of  the  target. 

Xt(t)  -  X  position  of  the  target  at  time  t. 

Yt(t)  -  Y  position  of  the  target  at  time  t. 

S„  -  X  component  of  target  speed. 

Si,  -  Y  component  of  target  speed. 

Rj  -  Steady-state  detection  range. 

X,„,Y„,X,(t),Y,(t),S„,S,,  -  Be  equivalent  variables  for  the  searcher. 

At  any  time  t,  the  X  and  Y  positions  of  the  target  and  searcher  can  be  written 
as  follows; 

X,(t)  =  X„  +  tSi.  and  Y,(t)  =  Y,,  +  tS„  (E.l) 

X.(t)  =  X„  +  tS„  and  Y.(t)  =  Y^  +  tS^  (E.2) 


2 


R(t)  =  ([X.(t)  -  XSt)f  +  [Y.(t)  -  YSt)fr 


(E.3) 


Combining  equations  (1),  (2),  (3)  and  setting  R(t)  =  R^  we  solve  for  t. 

t  =  -B  ±  (S'  -  4AC)‘’* 


where, 

A  =  -  2S„S„  +  S^,  +  S%  -  2S,^^ 

B  =  [2X«S„  +  2X,,S,.  -  2X^8..  -  2X,„S.J  +  [2Y„S^  +  2Y,„S,,  -  2Y„S,,  -  2Y„S^] 

C  =  +  r,,  -  2XJC,„  +  +  Y^,„  -  2Y„Y,.  -  Rd* 

The  time  to  steady-state  detection  is  the  minimum  of  the  positive  values  of 
t  obtained  from  equation  (E.4). 

Four  special  solution  cases  exist  in  which  the  time  to  steady-state  detection 
must  be  investigated  separately.  They  are:  when  both  values  of  t  are  negative, 
when  one  t  value  is  positive  and  the  other  is  negative,  when  (B*  -  4AC)  <  0,  and 
when  A=0. 

In  the  case  where  both  t  values  are  negative,  the  steady-state  detection 
opportunity  has  already  passed  and  detection  does  not  occur. 

When  one  t  value  is  positive  and  the  other  is  negative,  the  range  between 
the  searcher’s  and  target’s  initial  positions  is  less  than  the  steady-state  detection 
range.  As  a  result,  detection  occurs  and  the  time  to  detection  is  0. 

In  the  two  cases  (B^  -  4AC)  <  0  and  A  =  0  the  range  between  the  searcher 
and  the  target  never  decreases  to  the  steady-state  detection  range.  As  a  result, 
detection  does  not  occur. 


APPENDIX  F 


FORTRAN  CODE  FOR  SUBTRAN 

The  fortran  code  for  SUBTRAN  is  divided  into  the  following  parts: 

Main  Program  (SUBTRAN) 

Input  Subroutine  (IN) 

Propagation  Loss  Subroutine  (PLOSS) 

Target  Course  Subroutine  (TGCR) 

Target  Transient  Subroutine  (TGTR) 

Searcher  Course  Subroutine  (SRCR) 

Transient  Course  Subroutine  (TRCR) 

Steady-State  Detection  Subroutine  (SSDET) 

Transient  Detection  Subroutine  (TRDET) 

Lambda-Sigma  Jump  Subroutine  (LSJ) 

Position  Subroutine  (POSIT) 

Boundary  Reflection  Subroutine  (REFLCT) 

Track  Subroutine  (TRACK) 

Output  Subroutine  (OUT) 

Hold  Screen  Subroutine  (HSCRN) 

Uniform  Random  Number  Generator  Subroutine  (LRNDPC) 

Normal  Random  Number  Generator  Subroutine  (LNORPC) 

Gamma  Random  Number  (ienerator  Subroutine  (LGAMPC) 


74 


oooooo  oon  non  ooooooo 


cccccccccccccccccccccccccccccc 
c  c 

C  PROGRAM  SUBTRAN  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

PROGRAM  SUBTRAN 

*  SUBTRAN  requires  the  subroutines:  FALTR,  HSCRN,  IN, 
LSJ,  OUT,  PLOSS,  POSIT,  REFLCT,  SRCR,  SSDET,  TGCR, 
TGTR,  TRACK,  TRCR,  TRDET,  LRNDPC,  LNORPC,  LGAMPC. 

Declare  variables  and  dimension  arrays. 


CHARACTER*!  ANS 
DOUBLE  PRECISION  DSEED 

INTEGER  N,  NPpss,  NPptr,  NPc,  NPtc,  Ftr,  Fftr,  Ftw, 

1  Nc,  Nd,  Nz,  Nm,  Fstc,  Fsc,  Fssd,  Ftd,  HFsc 

REAL  L,  Tmax,  RNpss<50),  Ppss(50),  RNptr(50), 

1  Pptr<50),  St,  Ss,  Cc<50),  Pc<50),  Ttc<50), 

1  Ptc(50),  RPtr,  RPftr,  RPlsJ ,  SDlsJ ,  FOMss,  FOMtr, 

1  TV,  Rd(7),  T,  Xs,  Ys,  Xt,  Yt,  U(2),  Ct,  TCt ,  TTt, 

1  Cs,  TCs,  Xf,  Yf,  TTf,  ElsJ ,  Tlsj ,  Rssd,  Cr,  Tr, 

1  Tmin,  Tssd,  Td(5000),  HXs,  HYs,  HCs,  HTCs,  HTh 

.Read  input  data  from  the  user-defined  input  data  file. 

CALL  INCDSEED,  L,  N,  Tmax,  NPpss,  RNpss,  Ppss,  NPptr, 

1  RNptr,  Pptr,  St,  Ss,  NPc,  Cc,  Pc,  NPtc,  Ttc, 

1  Ptc,  Ftr,  RPtr,  Fftr,  RPftr,  RPlsJ,  SDlsJ, 

1  FOMss,  FOMtr,  Ftw,  TV) 

.Ask  if  input  parameters  are  correct. 

VRITEC*,  14) 

READ(*, 15) ANS 

IF((ANS  .NE.  ’Y’)  .AND.  <ANS  .NE.  ’y’))THEN 
VRITE(*, 24) 

CALL  HSCRN 
GOTO  60 
END  IF 
VRITE(*, 34) 


If  the  track  width  is  not  specified  by  the  user 
(Ftw=0)  then  compute  the  searcher’s  mean  track  width 
<TV)  in  nautical  miles  using  the  steady  state  figure 
of  merit  (FOMss)  and  the  propagation  loss  curve. 

IF(Ftw  .EQ.  0)THEN 

CALL  PLOSS (NPpss,  RNpss,  Ppss,  FOMss,  Rd) 
TV=2*Rd(l) 

END  IF 
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c 

C . start  counter  for  number  of  simulation  replications 

C  (Nc).  Initialize  detection  counter  (Nd),  detection  at 

C  time  zero  counter  (Nz),  and  detection  not  before  Tmax 

C  counter  (Nm)  . 


C 

Nc=l 

Nd=0 

Nz=0 

Nm=0 

C 


C . Assign  the  initial  value  of  the  simulation  time  <T)  in 

C  hours  to  zero.  This  is  the  starting  place  for 

C  successive  runs.  Set  the  searcher's  transient  course 

C  flag  (Fstc)  to  zero. 

C 

10  T=0 

Fstc=0 

C 

C . Set  the  initial  position  of  the  searcher  (Xs,Ys)  in 

C  the  lower  left-hand  corner  of  the  search  area  one-half 

C  a  track  width' s  distance  from  the  edges  of  the  search 

C  area. 

C 

Xs=TV/2 

Ys=TV/2 


Set  the  initial  position  of  the  target  (Xt.Yt) 
uniformly  distributed  In  the  search  area. 

CALL  LRNDPCCDSEED,  U.  2) 


Xt=L*U(l) 

Yt=L*U<2) 

C 

C . Set  initial  target  parameters  of  course  in  degrees 

C  (Ct),  time  to  next  course  change  in  hours  (TCt),  time 

C  to  next  transient  in  hours  (TTt).  Target’s  speed  is 

C  assumed  constant.  If  the  target’s  transient  flag  (Ftr) 

C  is  1,  transients  are  allowed  in  the  simulation. 


C 

CALL  TGCR<DSEED,  NPc,  Cc,  Pc,  NPtc,  If''.  Ptc,  Ct ,  TCt) 
IF(Ftr  .EQ.  DTHEN 

CALL  TGTRCDSEED,  RPtr,  TTt) 

ELSE 

TTt=2*Traax 
END  IF 
C 


C . Set  the  searcher’s  course  flag  (Fsc)  to  1. 

C  Set  initial  searcher  parameters  of  course  in  degrees 

C  <Cs),  time  to  next  course  change  n  hours  (TCs) . 

C  Searcher’s  speed  is  assumed  constant. 

C 
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Fsc=l 

CALL  SRCR(Fsc,  Xs,  Ys,  Ss,  TV,  L,  Cs,  TCs) 

C 

C . If  false  transient  detections  are  allowed  (Fftr=l), 

C  then  set  the  Initial  false  transient  parameters  of 

C  bearing  to  the  transient  in  degrees  CCf),  time  to  next 

C  false  transient  (TTf). 

C 

IFCFftr  .EQ.  DTHEN 

CALL  FALTR<DSEED,  RPftr,  L  Xf,  Yf,  TTf) 

ELSE 

TTf=2*Tmax 
END  IF 

. Set  initial  signal  error  level  (ElsJ )  and  time  time  to 

next  error  level  fluctuation  (TlsJ )  . 

CALL  LSJ (DSEED,  RPlsj .  SDlsj ,  Elsj ,  TlsJ ) 

. Compute  the  steady  state  detection  range. 

FOM=FOMss+ElsJ 

CALL  PLOSSCNPpss,  RNpss,  Ppss,  FOM,  Rd) 

Rssd=Rd ( 1 ) 

. Compute  the  time  In  hours  until  the  target  must  be 

reflected  off  the  search  area  boundry  (Tr)  If  present 
course  and  speed  are  maintained.  Also  compute  the 
course  after  reflection  (Cr). 

0  CALL  REFLCT(L,  Ct ,  St,  Xt,  Yt,  Cr,  Tr) 

. Compute  the  minimum  time  (Train)  until  the  next  event. 

30  Tmin=MIN( (Tmax-T) , TCt, TTt, TCs,  TTf, TlsJ , Tr) 

C 

C . Check  to  see  If  steady  state  detection  occurs  before 

C  Train.  If  so  stop  simulation  run  {Fssd=l).  Time  to 

C  steady  state  detection  (Tssd)  Is  computed  If  detection 

C  occurs. 

C 


CALL  SSDETCXt, 

Yt,  St, 

Ct,  Xs, 

Ys, 

Ss,  Cs,  Rssd,  Fssd, 

1 

Tssd) 

IF  < (Fssd 

.  EQ. 

1)  .AND. 

(Tssd 

.  LE. 

Train) ) THEN 

IF(  (T 

.  EQ. 

0)  .AND. 

(Tssd 

.  EQ. 

0)  )THEN 

Nz  = 

Nz  +  1 

GOTO  50 

ELSE 

GOTO  40 
END  IF 
END  IF 
C 
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C . Calculate  the  new  searcher  and  target  positions  based 

C  on  Tmln. 

C 

CALL  POSITCTmln,  Xt ,  Yt,  Ct,  St,  Xs,  Ys,  Cs,  Ss> 

C 

C . Decrement  the  time  values  and  check  to  see  what  the 

C  next  event  Is  . 


T=T+Tmln 

TCt=TCt-Tmin 

TTt=TTt-Tmln 

TCs=TCs-Tmln 

TTf=TTf-Tmln 

TlsJ=Tlsj-Tmln 

Tr=Ti: — Tmln 


C . Route  to  various  updates  based  on  which  of  the 

C  following  events  Is  due  to  occur. 

C  1.  Maximum  search  time  Is  up.  [T  >=  Tmaxl 

C  2.  Target  course  change.  [ TCt  =  03 

C  3.  Target  transient.  [ TTt  =  03 

C  4.  Searcher  course  change.  t TCs  =  03 

C  5.  False  transient.  [ TTf  =  03 

C  6.  Signal  error  fluctuation.  CTlsj  =  03 

C  7.  Target  boundry  reflection.  [Tr  =  03 

C 

C  —  Update  1. 


IF(T  .GE.  Traax)THEN 
Nm=Nm+ 1 
GOTO  50 
END  IF 
C 

C  —  Update  2. 

C 

IFCTCt  .EQ.  0)THEN 

CALL  TGCRCDSEED,  NPc,  Cc,  Pc,  NPtc,  Ttc,  Ptc,  Ct, 
1  TCt) 

GOTO  20 
END  IF 
C 

C  —  Update  3. 

C 

lF<TTt  . EQ.  0)THEN 

CALL  TGTRCDSEED,  RPtr,  TTt) 

FOM=FOMtr+ElsJ 

CALL  PLOSSCNPptr,  RHptr,  Pptr,  FOM,  Rd) 

CALL  TRDET(Xt,  Yt,  Xs.  Ys,  Rd,  Ftd) 

IFCFtd  .NE.  DTHEN 
GOTO  30 
ELSE 

Fstc= 1 
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CALL  TRCR(Fsc,  Fstc,  Xt,  Yt ,  Xs,  Ys,  Ss,  L,  TV, 
HFsc,  HXs,  HYs,  HCs,  HTCs,  HTh,  Cs, 
TCs) 

GOTO  30 
END  IF 
END  IF 

—  Update  4. 

IFCTCs  .EQ.  0)THEN 

IF(Fstc  .EQ.  0)THEN 

CALL  SRCRCFsc,  Xs,  Ys,  Ss,  TV,  L,  Cs,  TCs) 

GOTO  30 
ELSE 

CALL  TRCRCFsc,  Fstc,  Xt,  Yt ,  Xs,  Ys,  Ss,  L,  TV, 
1  HFsc,  HXs,  HYs,  HCs,  HTCs,  HTh,  Cs, 

1  TCs) 

GOTO  30 
END  IF 
END  IF 

—  Update  5. 

IFCTTf  .EQ.  0)THEN 
FOM=FOMtr+Elsj 

CALL  FLOSS (NPptr,  RNptr,  Pptr,  FOM,  Rd) 

CALL  TRDET<Xf,  Yf,  Xs,  Ys,  Rd,  Ftd) 

IF<Ftd  .NE.  DTHEN 

CALL  FALTRCDSEED,  RPftr,  L,  Xf,  Yf,  TTf) 

GOTO  30 
ELSE 

Fstc=l 

CALL  TRCR(Fsc,  Fstc.  Xf,  Yf,  Xs,  Ys,  Ss,  L,  TV, 
1  HFsc,  HXs,  HYs,  HCs,  HTCs,  HTh,  Cs, 

1  TCs) 

CALL  FALTRCDSEED,  RPftr,  L,  Xf ,  Yf ,  TTf) 

GOTO  30 
END  IF 
END  IF 

Update  6. 

IF(TlsJ  . EQ.  0)THEN 

CALL  LSJ (DSEED,  RPlsJ ,  SDlsJ ,  Elsj ,  Tlsj ) 
FOM=FOMss+ElsJ 

CALL  PLOSS<NPpss,  RNpss,  Ppss,  FOM,  Rd) 

Rssd=Rd  ( 1 ) 

GOTO  30 
END  IF 

—  Update  7. 
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IFCTr  . EQ.  0)THEN 
Ct=Cr 
GOTO  20 
END  IF 
C 

C . Store  times  to  detection  (Td) .  Update  detection 

C  counter  (Nd) .  Update  replication  counter  (Nc).  Route 

C  for  next  replication  run  If  required.  Print  to  the 

C  screen  a  message  every  100  replications. 

C 

40  Nd=Nd+l 

Td(Nd)=T+Tssd 
50  Nc=Nc+l 

lF<{INT<<Nc-l)/<N/10)>*(N/10>)  . EQ.  (Nc-1>)THEN 
VRITE(*, 44) (Nc-1) 

END  IF 

IF(Nc  . LE.  N)GOTO  10 
C 

C . Compute  the  output  statistics  and  print  to  the  user's 

C  terminal  and  output  file. 

C 

CALL  0UT(Nd,  Nz,  Nra,  Ftw,  TV,  Tmax,  Td) 

C 

C . Print  a  closing  message  to  the  user’s  terminal. 

C 

VRITE(*, 54) 

C 

C .  Read  format  statements. 

C 

15  F0RKAT<A) 

C 

C .  Write  format  statements. 

C 

14  FORMAT (23 (/),  IX,  ' Are  the  input  parameters  correct?’, 

1  ’ (Y/N) ’ /) 

24  FORMAT <23 </), 27X, ’ EXECUTION  TERMINATED  BY  USER’/ 

1  27X,28(’-’  ),  13(/)) 

34  FORMAT (23 (/) , IX, ' Executing. ...’/) 

44  FORMAT (5X, 14, ’  replications  complete.’) 

54  FORMAT(23  (/)  ,  32X,  6  (’ *  •)/32X,’)»:  SUBTRAN  *’ / 

1  32X,  8  (’  *  ’  )/// 

1  30X,  ’  Simulat ion  Complete’  , 9 (/) ) 

60  END 
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cccccccccccccccccccccccccccc 

c 

C  SUBROUTINE  IN 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 


c  c 
c 
c 
c 

c  c 


c 

c 

c 

c 

c 

c 

c 

c. 

c 


SUBROUTINE  INCDSEED,  L,  N.  Tmax,  NPpss,  RNpss,  Ppss, 

1  NPptr,  RNptr,  Pptr,  St,  Ss,  NPc,  Cc,  Pc, 

1  NPtc,  Ttc,  Ptc,  Ftr,  RPtr,  Fftr,  RPftr, 

1  RPlsJ,  SDlsj,  FOMss,  FOMtr,  Ftw,  TV) 

Variables  passed  to  INPUT  are;  none 

Variables  returned  by  INPUT  are:  all  in  subroutines 

call  list 

*  IN  requires  the  subroutine  HSCRN. 

.Define  variables.  Dimension  arrays. 


CHARACTER*12  FILEIN,  FILEOUT 
DOUBLE  PRECISION  DSEED 

INTEGER  N,  NPpss,  NPptr,  NPc,  NPtc,  Ftr,  Fftr,  Ftw, 

1  I  HOLD 

REAL  L,  Tmax,  RNpss (50),  Ppss (50),  RNptr (50), 

1  Pptr(50),  St,  Ss,  Cc(50),  Pc(50),  Ttc(50), 

1  Ptc(50),  RPtr,  RPftr,  RPlsj ,  SDlsJ ,  FOMss,  FOMtr, 

1  TV 


Write  to  the  user's  terminal  an  Intoructlon  message. 

WRITE(*, 14) 

CALL  HSCRN 


. Read  the  name  of  the  input  data  file  from  the  user's 

ermlnal.  Assign 

the  name  of  the  output  file  to  be  used  to  store  the 
nput  data. 

VRITE(*, 24) 

READ(*, 15)F1LEIN 
F I LEOUT= ' SUBTRAN .  OUT ' 


Open  the  user-defined  input  data  file  and  the  assigned 
output  data  file. 

OPENdO,  FILE=FILEIN) 

OPEN(20, FILE=FILEOUT) 

C 

C . Read  the  15  data  fields  from  the  input  data  file. 

C 

C  —  Data  field  1.  I  DSEED, L, N, TmaxJ 
C 


m 


81 


c 

c  — 
c 

c 

c  — 
c 


10 

c 

c  — 
c 

c 

c  — 
c 


20 

c 

c  — 
c 

c 

c  — 
c 

c 

c  — 
c 


30 

c 

c  — 
c 

c 

c  — 
c 


40 

c 

c  — 
c 

c 


READCIO, 25)DSEED, L, N, Tmax 
Data  field  2.  L NPpssl 
READdO,  35)NPpss 

Data  field  3.  E RNpss (NPp) , Ppss (NPp) ] 

DO  10  I=l,NPpss 

READ (10,45) RNpss  < I ) , Ppss  < I ) 
CONTINUE 

Data  field  4.  [ NPptrl 

READdO,  35)NPptr 

Data  field  5.  f RNptr (NPp) , Pptr ( WPp) ] 

DO  20  1=1 , NPptr 

READ (10,45) RNptr ( I ) , Pptr ( I ) 
CONTINUE 

Data  field  6.  tSt,Ss] 

READdO,  55)St,  Ss 
Data  field  7.  t NPcl 
READdO,  35)  NPc 

Data  field  8.  [ Cc (NPc) , Pc (NPc) I 

DO  30  1=1, NPc 

READ (1 0 . 45 ) Cc  (I ) , Pc  (  I ) 

CONTINUE 

Data  field  9.  [ NPtcI 

READdO,  35)NPtc 

Data  field  10.  t Ttc (NPtc) , Ptc (NPtc ) ] 

DO  40  1=1, NPtc 

READdO,  45)Ttc(I),  Ptc  ( I  ) 

CONTINUE 

Data  field  11.  [Ftr,RPtr] 

READdO,  65)Ftr,  RPtr 
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I 


—  Data  field  12.  tFftr.RPftrl 


READdO,  65)Ff  tr,  RPftr 


—  Data  field  13.  [  RPlsJ , SDlsJ ] 


READdO,  55)RPlsJ  ,  SDlsJ 


—  Data  field  14.  [  FOMss, FOMtr] 


READ (10,55) FOMss ,  FOMtr 


—  Data  field  15.  [  Ftw, TV] 


READdO,  65)Ftw,  TV 


Print  to  the  user's  terminal  and  the  output  data  file 
a  verification  of  the  input.  There  are  three  screens 
of  information. 

1.  Miscellaneous  data. 

2.  Propagation  loss  curve  -  steady  state 

(points) . 

3.  Propagation  loss  curve  -  transient 

(points) . 

4.  Target’s  course  and  time  to  course 
change  curves  (points) . 


Screen  1. 


VRITEO*:,  4) 

VRITE(*, 34) 

VRITE(20, 34) 

VRITE(*, 44)N. (L**2) , DSEED, Tmax 
VRITE(20, 44)N, (L**2 ) , DSEED, Tmax 
:F(Ftw  .EQ.  DTHEN 
VRITE(*, 54)TV 
VRITE(20, 54)TV 
ELSE 

VRITE(*, 64) 

VRITE(20, 64) 

END  IF 

VRITE(*, 74)RPlsJ , SDlsJ , Ftr, RPtr, Fftr, Ss, RPftr, St 
VR ITE (20, 74) RPlsJ , SDlsJ , Ftr , RPtr , Ff tr , Ss, RPf  tr , St 
CALL  HSCRN 


Screen  2, 


VRITE(*, 4) 
VRITE(*, 84)FOMss 
VRITE(20, 84) FOMss 
IHOLD=0 

DO  50  1=1, NPpss 


IHDLD=IH0LD+1 

WRITE <t, 94)RNpss(l ) , Ppss ( I ) 

VR I TE (20 , 94 ) RNpss ( I ) , Ppss ( I ) 

IF((]  .EQ.  15)  .OR.  (1  .EQ.  30)  .OR.  (1  .  EQ.  45)) 

1  THEN 

VRITEd;,  *) 

WRITE (K, *) 

CALL  HSCRN 
WRITE (^, 104) 

IHOLD=0 
END  IF 

50  CONTINUE 

IF<<I  .NE.  16)  .OR.  (I  .NE.  31)  .OR.  (I  .  NE.  46))THEN 

DO  60  1=1, 19- I HOLD 
WRITE(*. *) 

60  CONTINUE 

CALL  HSCRN 
END  IF 

WRITE<20. 154) 

C 

C  —  Screen  3. 

C 

WRITE(1:,  4) 

WRITECt:,  86)F0Mtr 
WRITE(20. 86)FOMtr 
IHOLD=0 

DO  70  1=1, NPptr 
IHOLD=IHOLD+l 

WRITEC*,  94)RNptr(I )  ,  Pptrd  ) 

WR I TE (20 , 94 ) RNptr ( I ) , Pptr  ( I ) 

IF((I  .EQ.  15)  .OR.  (I  .EQ.  30)  .OR.  (I  .  EQ.  45)) 

1  THEN 

WRITE(*,  *) 

WRITE(*,  *) 

CALL  HSCRN 
WRITE(*, 104) 

IHOLD=0 
END  IF 

70  CONTINUE 

IF((I  .NE.  16)  .OR.  (I  .NE.  31)  .OR.  (I  .NE.  46))THEN 
DO  80  1=1, 19-IHOLD 
WRITE(*, *) 

80  CONTINUE 

CALL  HSCRN 
END  IF 

WRITE(20, 154) 

C 

C  --  Screen  4. 

C 

WRITE(*, 4) 

WRITE (*, 114) 
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WRITE (20, 114) 

IHOLD=0 

DO  90  1  =  1,  MAXCNPc, NPtc) 

IHOLD=IHOLD+l 

IF<<I  .LE.  NPc)  .AND.  (I  .LE.  NPtc))THEN 
VRITEC^;,  124)Cc(I),Pc(I),Ttc(I).Ptc(I) 
VRITE(20,  124>Cc(I)  ,  Pc(I>  ,  Ttc(I  )  ,  PtcCI) 

GOTO  100 
END  IF 

IF(I  .GT.  NPtc)THEN 

WRITE (*, 134)Cc(I),Pc(I) 

WRITE <20, 134)Cc<l),Pc(I> 

GOTO  100 
END  IF 

IF (I  .GT.  NPc) THEN 

WRITE <*, 144)Ttc<l) , Ptc ( I ) 

WRITE <20,  144)Ttc  <I ) , Ptc  < I ) 

GOTO  100 
END  IF 

100  1F<<I  .EQ.  18)  .OR.  <I  .EQ.  36))THEN 

WRITEOt:,  ♦) 

CALL  HSCRN 
WRITE  <*,  104) 

IHOLD=0 
END  IF 

90  CONTINUE 

1F<<I  .NE.  19)  .OR.  <I  .NE.  37))THEN 
DO  110  I=1,20-IHOLD 
WRITE  C*;,  *) 

110  CONTINUE 

CALL  HSCRN 

END  IF 

WRITE<20, 154) 

C 

C .  Read  format  statements. 

C 

15  FORMAT <A) 

25  FORMAT  <23  </ ) , 2X , D13 . 0 , 2X,  F13 . 0 ,  2X ,  1 13 , 2X, F13 . 0 ) 

35  FORMAT<7 </) , 2X,  1 13, 8 </) ) 

45  FORMAT(2X, F13.  0,  2X, Fi3.  0) 

55  FORMAT <8 </) ,2X,F13. 0,2X,F13.  0) 

65  FORMAT <8 </) , 2X,  1 13, 2X, F13. 0) 

C 

C .  Write  format  statements. 

C 

4  FORMAT <23 </) ) 

14  format <23 </),  20X,  '  Submarine  Transient  Detection  ', 

1  ’ Simulation’ , 5 </), 32X,  8 <•  *  ’)/32X, 

1  SUBTRAN  *’/32X,8<’^  ’),5</),29X, 

1  ’J.  Brad  Kratovll  LT, USN’ /35X, ’ 26  July  ’, 

1  ' 1988' , 6 (/) ) 

24  FORMAT <23 </), IX, ’ ENTER  THE  NAME  OF  THE  INPUT  DATA  ’, 
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1  ’ FILE. ' /> 

34  F0RMAT(2OX.  IOC  ^  '  )/29X,  ’  *  SUBTRAN  INPUT  *’/ 

1  29X,10C:*:  ’)> 

44  FORMAT (//T9, ■ SIMULATION  PARAMETERS’ , T50, ' SEARCH 
1  'PARAMETERS’/  T9, 21  ( ' - ’  ) , T50 ,  17  ( ’ - ’  ) / 

1  T9,  ’  Repl Icat ions  =  ' ,  1 6 , T50 ,  ’ Search  Area  =  ’, 

1  F9. 2/T9, ’ Random  Seed  =  ’,F9.2,T50, 

1  ’Search  Time  =  ’,F9.2> 

54  FORMAT (T50, ’ Track  Width  =  •,F9.2) 

64  FORMAT (T50, ’ Track  Width  =  calc.’) 

74  FORMAT (T9, ’ TRANSIENT  PARAMETERS’ /T9 , 20 (’-’), T50 , 

1  ’SIGNAL  ERROR  PARAMETERS’ /TIO, 

1  ’  ( 1  -  Yes  /  0  -  No) ’  . T50 , 24  < ’ - ’  ) /T50 , 

1  'Rate  Par.  =  ’ , F5. 2/T9, ’ Real  Transients :’, T50 , 

1  'Std.  Dev.  =  ’, F5. 2/TlO, ’ Allowed  =  ’,12/TlO, 

1  'Rate  Par.  =  ’ , F7 . 4 . T50 , ’ SPEED  PARAMETERS’ /T9 , 

1  'False  Transients; *, T50 , 16 <’-’) /TIO , 

1  'Allowed  =  ’,I2,T50,’S.  Speed  =  ’,F5.2/T10, 

1  'Rate  Par.  =  ’ , F7 . 4 , T50 , ’ T .  Speed  =  ’,F5.2, 

1  4</)) 

84  FORMAT (12X, ’ STEADY-STATE  FIGURE  OF  MERIT  /  ’, 

1  'PROPAGATION  LOSS  PARAMETERS’ /12X, 

1  59 ('-’) //22X, ’ Steady-state  Figure  of  ', 

1  ’ Merit  =  ’ , F6. 2/) 

86  FORMAT ( 13X, ’ TRANSIENT  FIGURE  OF  MERIT  /  PROPAGATION  ’ 
1  'LOSS  ’,’ PARAMETERS’ /13X, 56 C -’ )//  23X, 

1  'Transient  Figure  of  ’,’ Merit  =  ’,F6.2/) 

94  FORMAT (T 13, ’ Range  Point  ;  ',F6.2,T50, 

1  'Prop.  Loss  :  ’,F6.2) 

104  FORMAT <8 (/)) 

114  FORMAT (28X, ’ TARGET  COURSE  PARAMETERS’ /28X. 24 (’-’) //) 
124  FORMAT (5X, ’ Course:  ' , F4. 0, 3X, ’ Prob:  ’,F4.2,T45, 

1  'Time  to  Change:  ’ , F6 . 2 , 3X, ’ Prob;  ’,F4.2> 

134  FORMAT (5X, ’ Course:  ' , F4 . 0 , 3X , ’ Prob :  ’,F4.2) 

144  FORMAT <T45, ’ Time  to  Change:  ’ , F6 . 2 , 3X, ’ Prob;  ’,F4.2) 
154  F0RMAT<3(/)) 

END 
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C  SUBROUTINE  FLOSS  C 

c  c 

cccccccccccccccccccccccccccccc 
c 

SUBROUTINE  FLOSS CNPp,  Rp.  Pp,  FOM,  Rd> 

Variables  passed  to  FLOSS  are;  NPp,  Rp(50),  Pp(50), 

FOM 

Variables  returned  by  FLOSS  are;  Rd(7> 

.  Declare  variables.  Dimension  arrays. 

INTEGER  NPp,  COUNT 

REAL  Rp<50),  Pp<50),  FOM,  Rd(7>,  SLOPE,  Rlnt 

. Initially  assign  the  detection  ranges,  Rd,  to  zero. 

DO  10  1=1,7 
Rd<I)=0 
0  CONTINUE 

. Compute  the  Intersection  of  the  propagation  loss  curve 

and  the  figure  of  merit. 

COUNT=l 

DO  20  1=1, NPp- 1 

SLOPE= (Pp ( I + 1 ) -Pp ( I ) ) / <Rp ( I + 1 ) -Rp ( I > ) 

Rlnt= ( CFOM-Pp  < I ) ) /SLOPE) +Rp ( I ) 

IF((Rlnt  .GE.  Rp(I))  .AND.  (Rlnt  . LE.  Rp ( I + 1 ) ) ) THEN 
Rd(COUNT)=Rlnt 
COUNT=COUNT+l 
IF (COUNT  .EQ.  8) GOTO  30 
END  IF 
20  CONTINUE 

30  END 
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ccccccccccccccccccccccccccccc 
c 

C  SUBROUTINE  TGCR 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  TGCRCDSEED,  NPc,  Cc,  Pc.  NPtc,  Ttc,  Ptc, 

1  Ct,  TCt) 

C 

C  Variables  passed  to  TGCR  are:  DSEED,  NPc,  Cc,  Pc, 

C  NPtc,  Ttc,  Ptc 

C  Variables  returned  by  TGCR  are:  Ct,  TCt 

C 

C  *  TGCR  requires  the  subroutine  LNORPC. 

C 

C .  Declare  variables.  Dimension  arrays. 

C 

DOUBLE  PRECISION  DSEED 
INTEGER  NPc,  NPtc 
REAL  Cc(50).  Pcc50),  Ttc<50>,  Ptc<50),  Ct,  TCt,  U<2>, 

1  SLOPE 

C 

C . Generate  two  Unlform<0,l)  random  variables  to  be  used 

C  In  extracting  the  new  course  and  time  to  course 

C  change. 

C 

CALL  LRNDPC (DSEED,  U,  2) 

C 

C . Extract  the  new  course  value  (Ct). 

C 

DO  10  1=1, NPc- 1 

SLOPE= (Pc(I+l)-Pc(I ) )/(Cc(I+l)-Cc(l ) ) 

IF((U(1)  .GE.  Pc(I))  .AND.  (U(l)  . LE.  Pc ( I + 1 ) ) ) THEN 
Ct= ( (U (l)-Pc (I ) )/SLOPE>+Cc ( I ) 

GOTO  20 
END  IF 

10  CONTINUE 

C 

C . Extract  the  new  time  to  course  change  (TCt). 

C 

20  DO  30  I=l,NPtc-l 

SLOPE= (Ptc ( I  + 1 ) -Ptc  < I ) ) / (Ttc ( I  + 1 ) -Ttc  ( I ) ) 

IF((U(2)  .GE.  Ptc(I))  .AND.  (U(2)  .  LE.  Ptc(I  +  l))) 

1  THEN 

TCt= ( (U (2 ) -Ptc ( I > ) /SLOPE) +Ttc ( I ) 

GOTO  40 
END  IF 

30  CONTINUE 

40  END 
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cccccccccccccccccccccccccccccc 
c  c 

C  SUBROUTINE  TGTR  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  TGTRCDSEED,  RPtr,  TTt) 

Variables  passed  to  TGTR  are:  DSEED,  RPtr 
Variables  returned  by  TGTR  are;  TTt 

*  TGTR  requires  the  subroutine  LGAMPC. 

Declare  variables. 

DOUBLE  PRECISION  DSEED 
REAL  RPtr, TTt , E 

Compute  the  time  to  next  transient  (TTt).  The  time 
between  transients  Is  exponentially  distributed  with 
mean  =  1/RPtr  in  hours. 

CALL  LGAMPC (DSEED, E, 1,  1.  ) 

TTt=E/RPtr 
END 


I 


I 
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cccccccccccccccccccccccccccccc 
c  c 

C  SUBROUTINE  SRCR  C 

c  c 

cccccccccccccccccccccccccccccc 
c 

SUBROUTINE  SRCR(C,  Xs,  Ys,  Ss,  TV,  L,  Cs,  TCs) 

Variables  passed  to  SRCR  are:  C,  Xs,  Ys,  Ss,  TV,  L 
Variables  returned  by  SRCR  are:  Cs,  TCs 

♦  The  variable  C  Is  generated  by  SRCR  when  It  Is  first 
run  and  Is  passed  back  to  SRCR  on  successive  runs. 

Declare  variables. 

INTEGER  C 

REAL  Xs,  Ys,  Ss,  TV,  L,  Cs,  TCs,  PI 
Assign  the  value  of  PI. 

PI=3. 1415926 

Set  up  the  program  routing  to  control  the  course. 


1F(  <C 

.  EQ. 

1) 

.OR. 

<C 

.EQ. 

5) ) GOTO 

10 

IF(  (C 

.  EQ. 

2) 

.  OR. 

<C 

.  EQ. 

6)  ) GOTO 

20 

IF<  (C 

.  EQ. 

3) 

.OR. 

<C 

.  EQ. 

7) )GOTO 

30 

IF(  (C 

.EQ.  4) 

.OR. 

(C 

.  EQ. 

8)  )G0TD 

40 

—  Searcher's  course  -  north. 

0  Cs=0 

TCs=  <L-TV)/Ss 
IF<C  .  EQ.  DTHEN 
C=2 
ELSE 
C=6 
END  IF 
GOTO  50 

—  Searcher’s  course  -  east  or  west  (following  north). 

20  IF(C  .EQ.  2)THEN 

Cs=90 

lF<Xs  .GE.  <L-<1.5*TV)))THEN 
TCs=  < (L- <TV/2) )-Xs)/Ss 
C=7 
ELSE 

TCs=TV/Ss 
C=3 
END  IF 
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ELSE 

Cs=270 

IF<Xs  .LE.  <1. 5*TV) )THEN 
TCs= (Xs- (TV/2 ) ) /Ss 
C=3 
ELSE 

TCs=TV/Ss 
C=7 
END  IF 
END  IF 
GOTO  50 

—  Searcher's  course  -  south. 

0  Cs=180 

TCs= (L-TV)/Ss 
IF(C  .EQ.  3)THEN 
C=4 
ELSE 
C=8 

END  IF 
GOTO  50 

Searcher’s  course  -  east  or  west  (following  south). 

40  IF(C  ,EQ.  4)THEN 

Cs=90 

IF(Xs  .GE.  (L-(1.5*TV)))THEN 
TCs= ( (L- (TV/2) )-Xs)/Ss 
C=5 
ELSE 

TCs=TV/Ss 
C=1 
END  IF 
ELSE 

Cs=270 

IF(Xs  .LE.  (1. 5*TV))THEN 
TCs= (Xs- (TV/2) )/Ss 
C=1 
ELSE 

TCs=TV/Ss 
C=5 
END  IF 
END  IF 
GOTO  50 
50  END 


91 


cccccccccccccccccccccccccccccc 
c  c 

C  SUBROUTINE  TRCR  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  TRCR<C,  Cl,  Xt,  Yt,  Xs,  Ys,  Ss,  L,  TV,  HC, 
1  HXs,  HYs,  HCs,  HTCs,  HTh,  Cs,  TCs) 

C 

C  Variables  passed  to  TRCR  are:  C,  Cl,  Xt,  Yt,  Xs,  Ys, 

C  Ss,  L,  TV 

C  Variables  returned  by  TRCR  are:  Cs,  TCs 

C 

C  *  The  variables  HC,  HXs,  HYs,  HCs,  HTCs,  HTh  are 

C  generated  by  TRCR  when  It  Is  first  run  and  are 

C  passed  back  to  TRCR  on  successive  runs. 

C 

C  *  TRCR  requires  the  subroutine  TRACK. 

C 

C . Declare  variables. 

C 

INTEGER  C,  Cl,  HC 

REAL  Xt,  Yt,  Xs,  Ys,  Ss,  L,  TV,  HXs,  HYs,  HCs,  HTCs, 

1  HTh,  Cs,  TCs 

C 


C . Set  up  routing  depending  on  Cl  value  as  follows: 


c 

1. 

Transient  first  detected. 

[ Cl=l] 

c 

2. 

Return  back 

along  the  transient  course. 

[ Cl=2] 

c 

3. 

Return  to  "on  hold" 

steady  state  search 

c 

position . 

t  Cl=3] 

c 

c 

4. 

Resume  steady  state 

search. 

I Cl=4] 

IF  (Cl 

.  EQ. 

DGOTO 

10 

IF(C1 

.  EQ. 

2) GOTO 

20 

IF(C1 

.  EQ. 

5) GOTO 

30 

p 

IF(C1 

.  EQ. 

4 ) GOTO 

40 

c  — 

p 

Route 

1. 

w 

10 

IF(C  . 

LE. 

8) THEN 

HXs=Xs 
HYs=Ys 
HCs=Cs 
HTCs=TCs 
HC=C 
C=9 
END  IF 

CALL  TJ?ACK  (Xt ,  Yt ,  Xs,  Ys,  Ss,  L,  Cs,  TCs) 
TCs=TCs- (TV/ (2*Ss) ) 

HTh=TCs 
Cl=2 
GOTO  50 
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Route  2. 


C 

C  " 

C 

20  IF(Cs  .LE.  180)THEN 
Cs=Cs+180 
ELSE 

Cs=Cs-180 
END  IF 
TCs=HTh 
Cl=3 
GOTO  50 

C 

C  —  Route  3. 

C 

30  CALL  TRACK (HXs,  HYs,  Xs.  Ys,  Ss,  L,  Cs,  TCs) 
TCs=SQRT( ( <HYs-Ys)**2)+< (HXs-Xs)**2> )/Ss 
Cl=4 
GOTO  50 

C 

C  —  Route  4. 

C 

40  C=HC 

Cs=HCs 
TCs=HTCs 
C1  =  0 

50  END 
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c  c 
c 
c 
c 

c  c 
c 


Variables  passed  to  SSDET  are:  Xt,  Yt,  St,  Ct,  Xs,  Ys, 

Ss,  Cs,  Rssd 

Variables  returned  by  SSDET  are:  Cl,  Tssd 
.  Declare  variables. 

INTEGER  Cl 

REAL  Xt,  Yt,  St,  Ct,  Xs,  Ys,  Ss,  Cs,  Rssd,  Tssd,  PI, 

1  HOLD,  Stx,  Sty,  Ssx,  Ssy,  A,  B,  C,  Ax,  Bx,  Cx, 

1  Ay,  By,  Cy,  Tssda,  Tssdb 

•Assign  the  value  of  PI. 

PI=3. 1415926 

.Check  to  see  If  the  searcher  is  already  within  steady 
state  detection  range. 

HOLD=SQRT  <  (  (Xs-Xt )  )  +  <  <Ys-Yt )  **2)  ) 

IF (HOLD  .LE.  Rssd) THEN 
Tssd=0 
Cl  =  l 
GOTO  20 
END  IF 

.Compute  the  target's  x  speed  component  (Stx)  and  y 
speed  component  (Sty) . 

Stx=St»SIN(Ct*PI/180) 

Sty=St*COS(CttP 1/180) 

.Compute  the  searcher's  x  speed  component  (Ssx)  and  y 
speed  component  (Ssy) . 

Ssx=Ss*SIN(Cs*PI/180) 

Ssy=Ss*C0S (Cs^P I /180 ) 

.Compute  the  constants  to  be  used  in  the  quadratic 
equation  for  solving  the  time  to  steady  state 
detection  (if  It  occurs). 


CCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

SUBROUTINE  SSDET  C 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  SSDET  (Xt,  Yt,  St,  Ct ,  Xs,  Ys,  Ss,  Cs,  Rssd, 
1  Cl,  Tssd) 


Ax= (Ssx**2) + (Stx**2)- (2*Ssx*Stx) 

Bx= (2tXstSsx)+ (2*Xt*Stx)- (2*XstStx)- (2*Xt*Ssx) 
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Cx= (Xs»*2)+ (XtK*2)- <2*Xs*Xt) 

Ay= (Ssy**2) + (Sty*^2 ) - <2*Ssy*Sty) 

By=  <2*Ys^Ssy )  +  (2*yt*Sty ) - (2*Ys*Sty ) - (2*Yt*Ssy > 

Cy= (Ys**2)  + <Yt**2)- (2*Ys*Yt) 

A=Ax+Ay 

B=Bx+By 

C=Cx+Cy- <Rssd**2) 

....Check  the  conditions  that  yield  a  solution  that  is  not 
real  or  finite. 

1.  A=0 

2.  (B**2)- (4*A*C)<0 
Condition  1. 

IF<A  .EQ.  0)GOTO  10 

—  Condition  2. 

HOLD=  <B**2)  -  (4:»:A*C) 

IFCHOLD  .LT.  0)GOTO  10 

.  .  .  .Compute  the  values  for  the  time  to  steady  state 
detection  (Tssd)  . 

Tssda=  <-B+SQRT (HOLD) ) / (2*A) 

Tssdb=  <-B-SQRT<HOLD) )/(2*A) 

....Check  for  negative  detection  times. 

1.  Both  Tssda  and  Tssdb  negative. 

2.  Tssda  positive,  Tssdb  negative. 

3.  Tssda  negative,  Tssdb  positive. 

4.  Both  Tssda  and  Tssdb  positive. 

—  Case  1. 

IF ((Tssda  . LT.  0)  .AND.  (Tssdb  . LT.  0)>GOTO  10 

—  Case  2. 

IF( (Tssda  .GE.  0)  .AND.  (Tssdb  .LT.  0))THEN 
Tssd=Tssda 
Cl  =  l 
END  IF 

Case  3. 

IF(  (Tssda  .LT.  0)  .AND.  (Tssdb  .GE.  0))THEN 
Tssd=Tssdb 
Cl  =  l 
END  IF 


95 


o  o 


Case  4 . 

lF(<Tssda  .GE.  0)  .AND.  <Tssdb  .GE. 
Tssd=MIN (Tssda, Tssdb) 

Cl  =  l 
END  IF 
GOTO  20 
10  C1=0 

20  END 


0) >THEN 


ooot-*  non  ooooooonon  onoo  noooon  oooooo 


ccccccccccccccccccccccccccccc 

c 

SUBROUTINE  TRDET  C 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  TRDETCXt,  Yt,  Xs,  Ys,  Rd,  C) 

Variables  passed  to  TRDET  are:  Xt,  Yt,  Xs,  Ys,  Rd(7) 
Variables  returned  by  TRDET  are:  C 

. Declare  variables.  Dimension  arrays. 

INTEGER  C 

REAL  Xt,  Yt,  Xs,  Ys,  Rd<7),  HOLD 

. Compute  the  current  range  between  the  searcher  and 

target  (HOLD) . 

HOLD=SQRT<  (  ( Xs-Xt )  )  +  (  (Ys-Yt>*t2>  ) 

. Check  to  see  if  the  current  range  is  within  transient 

detection  range.  The  code  (C)  for  a  detection  Is  1, 
and  the  code  for  no  detection  Is  0.  Two  cases  exist 
that  require  checking: 

1.  Direct  path  detection. 

2.  Convergence  zone  detection  (3  max). 

Case  1. 

IFCHOLD  .LE.  Rd(l>>THEN 
C=1 

GOTO  20 
END  IF 

Case  2. 

DO  10  1=2, 6, 2 

IF((HOLD  .GE.  Rd(l))  .AND.  (HOLD  . LE.  Rd(l+1)))THEN 
C=1 

GOTO  20 
END  IF 

0  CONTINUE 

. Indicate  that  no  detection  occurred. 

C=0 
20  END 
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ccccccccccccccccccccccccccc 
c 

C  SUBROUTINE  L3J 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  LSJ (DSEED,  RPlsj ,  SDlsJ ,  ElsJ ,  TlsJ ) 

Variables  passed  to  LSJ  are:  DSEED, RPlsJ , SDlsJ 
Variables  returned  by  LSJ  are:  ElsJ, TlsJ 

♦  LSJ  requires  the  subroutines:  LNORPC,  LGAMPC. 

Declare  variables. 

DOUBLE  PRECISION  DSEED 
REAL  RPlsj,  SDlsJ,  ElsJ,  TlsJ,  N,  E 

Calculate  the  Lambda  Sigma  Jump  error  fluctuation 
(ElsJ).  The  error  fluctuation  Is  normally  distributed 
with  mean  zero  and  standard  deviation  =  SDlsJ. 

CALL  LNORPC (DSEED,  N,  1) 


ElsJ=SDlsJ  :t:N 
C 

C . Calculate  the  time  to  the  fluctuation  (TlsJ).  The  time 

C  between  flucuations  is  exponentially  distributed  with 

C  mean  =  1/RPlsj . 


C 

CALL  LGAMPC (DSEED.  E.  1,  1.) 

TlsJ=E/RPlsJ 

END 


C  C  C 
C 
C 
C 

C  C  C 
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ccccccccccccccccccccccccccc 
c 

C  SUBROUTINE  FALTR 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  FALTR (DSEED,  RPftr,  L,  Xf ,  Yf ,  TTf) 

Variables  passed  to  FALTR  are:  DSEED,  RPftr,  L 
Variables  returned  by  FALTR  are:  Xf,  Yf,  TTf 

*  FALTR  requires  the  subroutines:  LRNDPC,  LGAMPC 

Declare  variables.  Dimension  arrays. 

DOUBLE  PRECISION  DSEED 
REAL  RPftr,  L,  Xf ,  Yf,  TTf,  U(2>,  E 

Compute  the  x  and  y  positions  of  the  false  transient 
(Xf,Yf).  The  position  Is  uniformly  distributed  over 
the  search  area. 

CALL  LRNDPC (DSEED,  U,  2) 

Xf=L)KU  (1) 

Yf=L*U<2) 

Compute  the  time  to  the  next  false  transient  (TTf). 
The  time  between  false  transients  Is  exponentially 
distributed  with  mean  =  1/RPftr  in  hours. 

CALL  LGAMPC (DSEED,  E,  1,  1.) 

TTf=E/RPftr 
END 


C  C  C 
C 
C 
C 

C  C  C 
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ccccccccccccccccccccccccccccc 

c 

SUBROUTINE  POSIT  C 

c 

ccccccccccccccccccccccccccccc 

SUBROUTINE  POSITCT,  Xt,  Yt,  Ct,  St,  Xs,  Ys,  Cs,  Ss) 

Variables  passed  to  POSIT  are:  T,  Xt,  Yt,  Ct,  St,  Xs, 

Ys,  Cs,  Ss 

Variables  returned  by  POSIT  are:  Xt,  Yt,  Xs,  Ys 
*  POSIT  updates  the  variables  Xt,  Yt,  Xs,  Ys 
. . . .Declare  variables. 

REAL  T,  Xt,  Yt,  Ct,  St,  Xs,  Ys,  Cs,  Ss,  PI 
....Assign  the  value  for  PI. 

PI=3. 1415926 

....Compute  target’s  x  position  (Xt)  and  y  position  (Yt) . 

Xt  =  Xt+  (T*St*S IN  (Ct:»:P  1/180)  ) 

Yt=Yt+ (T*St*COS(Ct*PI/180) ) 

....Compute  searcher’s  x  position  (Xs)  and  y  position 
(Ys)  . 

Xs=Xs+ (T*Ss*SIN(Cs*PI/180) ) 

Ys=Ys+  (TtSs:*COS(Cs*:P  1/180)  ) 

END 
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cccccccccccccccccccccccccccc 
c 

C  SUBROUTINE  REFLCT 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  REFLCT (L.  Ct,  St.  Xt,  Yt ,  Cr,  Tr) 

Variables  passed  to  REFLCT  are:  L,  Ct,  St,  Xt ,  Yt 
Variables  returned  by  REFLCT  are:  Cr,  Tr 

Declare  variables. 

REAL  L,  Ct,  St,  Xt,  Yt,  Cr,  Tr,  PI,  Stx,  Sty,  Ty,  Tx 
Assign  the  value  of  PI. 

PI=3. 1415926 

Compute  the  x  component  of  target  speed  (Stx) .  Compute 
the  y  component  of  target  speed  (Sty) . 


Stx=St*SIN(Ct*PI/180) 

Sty=St*COS(Ct:*:P  1/180) 

C  ; 

C . Compute  the  time  to  reflection  (Tr)  and  the  reflected 

C  course  (Cr)  for  the  five  cases: 

C  1.  Special  cases  -  Ct=0, 90, 180, 270, 360 

C  2.  0<Ct<90 

C  3.  90<Ct<180 

C  4.  180<Ct<270 

C  5.  270<Ct<360 

C 

C  —  Special  cases. 


C 

IF((Ct  . EQ.  0)  .OR.  (Ct  .EQ.  360))THEN 
Tr= (L-Yt)/Sty 
Cr=180 
GOTO  10 

END  IF 

IF(Ct  .EQ.  90)THEN 
Tr= (L-Xt)/Stx 
Cr=270 
GOTO  10 

END  IF 

IF(Ct  .EQ.  180)THEN 
Tr=- Yt/Sty 
Cr=0 
GOTO  10 

END  IF 

IF(Ct  .EQ.  270)THEN 
Tr=-Xt/Stx 


C  C 
C 
C 
C 

C  C 
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Cr=90 
GOTO  10 

END  IF 

—  Case  1 . 

lF((Ct  .GT.  0)  .AND.  <Ct  .LT. 
Ty=  <L-Yt)/Sty 
Tx=  <L-Xt)/Stx 
Tr=MlN<Ty,Tx) 

IF<Tr  .EQ.  Ty)THEN 
Cr=180-Ct 
ELSE 

Cr=360-Ct 
END  IF 
GOTO  10 

END  IF 

Case  2 . 

IF(<Ct  .GT.  90)  .AND.  <Ct  .  LT. 
Ty=-Yt/Sty 
Tx=  <L-Xt)/Stx 
Tr=MINCTy, Tx) 

IF<Tr  . EQ.  Ty)THEN 
Cr=180-Ct 
ELSE 

Cr=360-Ct 
END  IF 
GOTO  10 

END  IF 

Case  3 . 

IFC(Ct  .GT.  180)  .AND.  (Ct  . LT 
Ty=-Yt/Sty 
Tx=-Xt/Stx 
Tr=MIN<Ty, Tx) 

IF(Tr  . EQ.  Ty)THEN 
Cr=540-Ct 
ELSE 

Cr=360-Ct 
END  IF 
GOTO  10 

END  IF 

Case  4 . 

IF(<Ct  .GT.  270)  .AND.  <Ct  . LT 
Ty= (L-Yt)/Sty 
Tx=-Xt/Stx 
Tr=MIN(Ty, Tx) 


90) ) THEN 


180) )THEN 


270) )THEN 


360) ) THEN 
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ccccccccccccccccccccccccccccc 

c 

SUBROUTINE  TRACK  C 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  TRACK <Xt,  Yt,  Xs,  Ys,  Ss,  L,  Cs,  TCs) 

Variables  passed  to  TRACK  are:  Xt,  Yt,  Xs,  Ys,  Ss,  L 
Variables  returned  by  TRACK  are:  Cs,  TCs 

♦  TRACK  requires  the  subroutine  REFLCT. 

. . . . Declare  variables. 

REAL  Xt,  Yt,  Xs,  Ys,  Ss,  L,  Cs,  TCs,  PI,  HOLD,  Cr,  Tr 
....Assign  the  value  for  PI. 

PI=3. 1415926 

.  .  .  .Compute  the  approach  track  course  (Cs)  for  the 
following  cases: 

1.  Special  case  where  the  y  positions  are  equal. 

2.  (Xt  .GE.  Xs)  and  (Yt  .GT.  Ys) 

3.  (Xt  .LE.  Xs)  and  (Yt  .LT.  Ys) 

4.  (Xt  .LE.  Xs)  and  (Yt  .GT.  Ys) 

5.  (Xt  .GE.  Xs)  and  (Yt  . LT.  Ys) 

—  Case  1 . 

IF(Ys  .EQ.  Yt)THEN 
IF(Xt  .GT.  Xs)THEN 
Cs=90 
ELSE 

Cs=270 
END  IF 
GOTO  10 
END  IF 

Compute  the  common  value  for  the  angle  of  the  track 
between  the  searcher  and  target  (HOLD)  to  be  used  In 
cases  2-5. 

HOLD= (ATAN( (Xt-Xs)/(Yt-Ys) ) >*180/PI 

—  Case  2. 

IF((Xt  .GE.  Xs)  .AND.  (Yt  .GT.  Ys))THEN 
Cs=HOLD 
GOTO  10 
END  IF 
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10 


Case  3. 


IF((Xt  .LE.  Xs) 
Cs=H0LD+180 
GOTO  10 

END  IF 

.  AND. 

(Yt 

.  LT. 

Ys) )THEN 

Case  4 . 

IF<(Xt  .LE.  Xs) 
Cs=HOLD+360 
GOTO  10 

END  IF 

.  AND. 

<yt 

.  GT. 

Ys) )THEN 

Case  5 . 

IF<(Xt  .GE.  Xs) 

.  AND. 

(Yt 

.  LT. 

Ys) )THEN 

Cs=HQLD+180 
END  IF 


Call  the  REFLCT  subroutine  to  get  the  time  to  next 
course  change  <TCs) . 

CALL  REFLCT<L,Cs,Ss,Xs,Ys,Cr.Tr) 

TCs=Tr 

END 
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c  c 
c 
c 
c 

c  c 
c 


Variables  passed  to  OUT  are:  Nd,  Nz,  Nm,  Ftw,  TV, 

Tnax,  Td 

Variables  returned  by  OUT  are:  none 

*  OUT  requires  the  subroutine  HSCRN 

.Declare  variables.  Dimension  arrays. 

CHARACTER*12  FILEl 

INTEGER  Nd.  Nz,  Nm,  Ftw,  ICOUNT 

REAL  TV,  Tmax,  Td<5000),  Tdor(5000>,  SUMd,  SUMm, 

1  SUMld,  SUMlm,  SUM2,  SUM3,  SUM4,  SUMS,  SUMS,  SUM7, 

1  XBAR,  SD,  SDx,  XBARd,  SBAR,  RPe,  RS,  SUR(5000>, 

1  PD(IO)  ,  PDTdO) 

.Designate  the  output  file.  FILEl  -  holds  the  times  to 
detection. 

FILE1=' TIME. OUT* 

.Open  FILEl  for  writing  to.  The  other  output  file 
(SUBTRAN.DAT)  Is  already  open  from  the  subroutine  IN 
and  Is  designated  as  device  20. 

OPEN (30 , F I LE=F I  LEI ) 


CCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

SUBROUTINE  OUT  C 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  OUT(Nd,  Nz,  Nm,  Ftw,  TV,  Tmax,  Td) 


Vrlte  to  the  user's  screen  a  message  indicating 
computing  output  statistics. 

VRITEC*, 14) 

Order  the  times  to  detection  from  smallest  to  largest 
(Tord) . 

Tdor (l)=Td(l) 

DO  10  1=1, Nd 

DO  20  J=I-1, 1,-1 

IF(Td(I)  .LT.  Tdor<J))THEN 
Tdor ( J  +  1 ) =Tdor (J ) 

Tdor (J)=Td<I) 

ELSE 

Tdor(J+l)=Td(I) 

GOTO  10 
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END  IF 

20  CONTINUE 

10  CONTINUE 

C 

C . Compute  the  mean  (XBAR>  and  standard  deviation  (SD)  of 

C  the  times  to  detection.  For  the  cases  where  detection 

C  did  not  occur  before  Tmax,  assign  a  value  of  Tmax  for 

C  the  time  to  detection.  Cases  where  the  time  to 

C  detection  Is  zero  are  not  Included.  Also  Include  the 

C  95%  confidence  bounds  on  the  estimate  of  the  mean 

C  <SDx). 

C 

DO  30  1=1 , Nm 

Tdor (Nd+ I ) =Tmax 
30  CONTINUE 


SUMd=0 

SUMm=0 

SUMld=0 

SUMlm=0 

DO  40  1=1, Nd 

SUMd=SUMd+Tdor ( I ) 
SUMld=SUMld+ (Tdor ( I >  **2) 

40  CONTINUE 

DO  50  1=1, Nm 

SUMm=SUMra+Tdar  <Nd+ 1 ) 
SUMlm=SUMlm+ (Tdor (Nd+ 1 ) **2) 
50  CONTINUE 

XBAR= (SUMd+SUMm) / (Nd+Nm) 


SD=SQRT( ( (SUKld+SUMlra)-( (Nd+Nm) *(XBAR**2) ))/(Nd+Nm-l> ) 
SDx=l. 645*SD/( (Nd+Nm)**.5) 

C 

C . Construct  a  empirical  survival  function  (SUR)  from  the 

C  ordered  time  to  detection  data.  Under  the  assumption 

C  that  the  data  Is  from  an  exponential  distribution, 

C  perform  a  linear  regression  on  the  natural  log  of  SUR 

C  with  respect  to  the  ordered  times  Tdor.  The  slope  of 

C  the  resultant  line  Is  an  estimate  for  the  rate 

C  parameter  of  the  exponential  distribution  (RPe) .  1/RPe 

C  is  another  estimate  for  the  mean  time  to  detection. 

C  This  estimate  is  independent  of  the  detections  that 

C  occurred  at  time  zero  thus  those  values  are  not 

C  included.  Also  computed  is  the  R-squared  value  of  the 

C  linear  regression.  This  value  is  an  Indication  of  how 

C  good  the  exponential  assumption  Is. 


C 

SUM2=0 
SUM3=0 
SUM4=0 
ICOUNT=0 
DO  60  1=1, Nd 

SUR ( I ) =1- ^REAL ( I ) /REAL (Nd+Nm) ) 
1F(SUR(I)  .LE.  0)THEN 
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IC0UNT=IC0UNT+1 
SUMd=SUMd-Tdor  <.  I ) 

SUMld=SUMld- (Tdar ( 1 ) ) 

GOTO  60 
END  IF 

SUM2=SUM2+AL0G (SUR ( I ) > 

SUM3=SUM3+ (Tdor ( I ) »ALOG (SUR ( I >  > ) 

SUM4=SUM4+ (ALOG (SUR ( 1 ) ) **2> 

60  CONTINUE 

DO  70  1=1, 10 
DO  80  J=l, Nd 

IF( (1-SUR(J ) )  .GE.  (REAL(I>/10))THEN 
PD(I )=1-SUR(J) 

PDT(I)=Tdor(J) 

GOTO  70 
END  IF 

80  CONTINUE 

70  CONTINUE 

XBARd=SUMd/ (Nd- I COUNT ) 

SBAR=SUM2/ (Nd- I COUNT) 

RPe=(SUM3- (XBARd*SUM2) )/ 

1  (SUMld- ( (Nd- ICOUNT) » (XBARd**2 >  > ) 

SUM5=0 

DO  90  1=1, Nd 

IF(SUR(I)  .LE.  0) GOTO  90 

SUM5=SUM5+ ( (Tdor ( I )-XBARd>  ^ (ALOG (SUR ( I ) ) -SBAR) > 

90  CONTINUE 

SU146=SUMld-  (  (Nd-ICOUNT)»(XBARd**2)  ) 

SUM7=SUM4-( (Nd-IC0UNT)*(SBAR**2>  > 

RS=  (SUM5**2 )  /  (SUM6<:SUM7) 

C 

C .  Write  the  output  to  the  screen  and  the  two  output  data 

C  files. 

C 

VRITE(*, 114) 

WRITE(*, 24) 

VRITE(20, 24) 

VRITE(*, 34)N2 
WRITE(20, 34)Nz 
VRITE(»,  44)Nni 
VRITE(20,  44)Nnj 
IF(Ftw  .EQ.  0)THEN 
WRITE(*,54)TV 
WRITE(20,54)TW 
ELSE 

VRITE(*, *) 

END  IF 

VRITE(^, 64)XBAR, SDx, SD, (-l/RPe> , RS 
VRITE(20,64)XBAR,SDx,SD, (-1/RPe) ,RS 
WRITE  (*,74) 

WRITE  (20,74) 

DO  100  1=1,5 
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VR I TE ( * , 84 ) PD ( I ) , PDT ( I ) , PD ( I +5 ) . PDT ( 1 +5 ) 

VR ITE (20 , 84 ) PD ( 1 ) , PDT ( 1 ) , PD ( I +5 ) , PDT ( I +5 ) 

100  CONTINUE 

VR1TE(*, *> 

CALL  HSCRN 
¥R1TE<*,94) 

DO  110  1=1, (Nd+Nm) 

VRITEOO,  *)Tdor(l) 

110  CONTINUE 

WRITEC*, 104) 

CALL  HSCRN 
C 

C . Write  format  statements. 

C 

14  FORMAT (23 (/), IX, ’ Computing  output  statistics....’/) 

24  FORMAT(28X, 12 (' *  ’)/28X, ’*  SUBTRAN  OUTPUT  DATA  *’/ 

1  28X, 12 (’ *  ’ )/) 

34  FORMAT (5X, • NO.  OF  REPS.  WITH  DETECTION  TIME  OF  ', 

1  •  ZERO  . ’  ,  14) 

44  FORMAT (5X, *  NO.  OF  REPS.  WITH  DETECTION  TIME  ’, 

1  ’GREATER  THAN  Traax . ’,14) 

54  FORMAT (5X, ’ CALCULATED  TRACK  WIDTH  VALUE  ', 

1  ’ . ’  ,  F7. 2) 

64  FORMAT (/26X,  ’ TIME  TO  DETECTION  STAT 1ST  ICS’ /26X, 

1  28 (’-’)/  26X,’Mean  =  ’,F7.2,’  +/-’, 

1  F5.2,’  (95%  Cl ) ’ /26X, ’ Std.  Dev.  =  ’,F7.2/ 

1  26X,  ’  1/Lambda  =  ’ , F7 . 2/26X, ’ R-Squared  =  ’, 

1  F7.2/) 

74  FORMAT (TIO, ' Prob.  of  Det . ’ , T27 , ’ By  time’,T47, 

1  ’Prob.  of  Det. ’ ,T64, ’ By  t Ime' /TIO , 13 ( ’ - ’ )  , 

1  T27, 7  (’ -’  ) , T47, 13(’ -’  ) , T64, 7  (’ -’  ) ) 

84  F0RMAT(T15, F3. 1, T28, F6. 1, T52, F3. 1, T65, F6. 1) 

94  FORMAT (23 (/), IX, ’ Writing  to  output  file....’/) 

104  FORMAT (4X, ’ Output  data  is  located  in  the  following  ’, 
1  ' files: ’ /7X. ’ 1.  SUBTRAN.DAT  -  input  ’, 

1  ’parameters,  time  to  detection  ’, 

1  ’ statistics. ’ /7X, ’ 2.  TIME. OUT  -  times  to  ’, 

1  'detection.'//) 

114  FORMAT (23 (/) ) 

END 
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non  non  nnnn 


ccccccccccccccccccccccccccccc 

c 

C  SUBROUTINE  HSCRN 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  HSCRN 

Write  to  the  user’s  terminal  a  "how  to  continue" 
message. 

WRITE(*, 14) 

Perform  a  generic  READ  statement  to  hold  the  screen. 
READ<*. *) 

Write  format  statement. 

14  FORMAT <26X,  *  -  Press  <ENTER>  to  continue  -') 

END 


110 


n  n  n  n  o 


noon 


cccccccccccccccccccccccccccccc 
c  c 

C  SUBROUTINE  LRNDPC  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  LRNDPC (DSEED,  U,  N> 

C 

INTEGER*4  N.  I 
REAL  U(N) 

DOUBLE  PRECISION  D31M1.  DSEED,  D31 

d31nil  =  2**31-1 
d31  =  2**31 

DATA  D31M1/2147483647.  DO/ 

DATA  D31/2147483648.  DO/ 


DO  5  I = 1 , N 

DSEED  =  DMOD < 16807. DO*DSEED. D31M1) 
U(I)  =  DSEED/D31 
CONTINUE 
END 
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o  o  o  o  o 


cccccccccccccccccccccccccccccc 
c  c 

C  SUBROUTINE  LNORPC 

C 

ccccccccccccccccccccccccccccc 
c 

SUBROUTINE  LNORPC < ISEED, A, N) 

This  program  will  generate  a  vector  of  normal  random 
variables  using  according  to  the  sine-cosine  method. 

INTEGER  N, I.IND 
REAL  A<N) , U (2) , S, V, Ul, XSTAR 
DOUBLE  PRECISION  ISEED.  PI 
DATA  PI/3. 14159265358979D0/ 

DATA  IND/1/ 

IND  =  1 
DO  100  1=1, N 
IND  =  -IND 

IF  (IND. GE. 0)  GOTO  20 
10  CALL  LRNDPC< ISEED,  U,  2) 

S  =SQRT<-2*ALOG(U(l) )  ) 

V  =  2*PItU(2) 

XSTAR  =  StCOS(V) 

A(I)  =  StSINCV) 

GOTO  100 

20  ACn  =  XSTAR 
100  CONTINUE 
RETURN 
END 
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o  o  n 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c  c 

C  SUBROUTINE  LGAMPC  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  LGAMPC (DSEED,  G,  N.  A) 

DOUBLE  PRECISION  DSEED 
DIMENSION  G(N),  UN<2) 

DATA  ASAVE  /- 1 . / . XLL/- 1 ,  / 

IF(A  . LE.  1. )  GO  TO  1500 

C  -  GAMMA  GENERATOR  A>1.  SCHMEIZER’S  ALGORITHM 

IF  (A  .EQ.  ASAVE)  GO  TO  50 

ASAVE  =  A 

X1=0 

X2=0 

F1  =  0 

F2=0 

X3  =  A  -  1. 

R3X=X3-X3»ALOG (X3) 

D  =  SQRT(X3) 

IF  (D  .  GE.  X3)  GO  TO  10 

X2  =  X3  -  D 

XI  =  X2*(l.-1./D) 

XLL  =  1.  -  X3/X1 

FI  =  EXP(X3*AL0G(X1)  -X1+R3X) 

F2  =  EXP(X3*ALOG<X2)  -  X2  +  R3X) 

10  X4  =  X3  +  D 

X5  =  X4*(l.  +  1. /D) 

XLR  =  1.  -  X3/X5 

F4  =  EXP(X3*ALOG<X4)  -  X4  +  R3X) 

F5  =  EXP<X33t:ALOG(X5)  -  X5  +  R3X) 

PI  =  F2*(X3-X2) 

P2  =  F4*(X4-X3)  +  PI 

P3  =  F1*<X2-X1)  +  P2 

P4  =  F5*{X5-X4)  +  P3 

P5  =  X3-X2  -PI  +  P4 

P6  =  X4-X3-P2+P1  +  P5 

P7  =  (F2-F1>*(X2-X1)*. 5  +  P6 

P8  =  (F4-F5)»(X5-X4)*. 5  +  P7 

P9  =  -Fl/XLL  +  P8 

PIO  =  F5/XLR  +  P9 

C  VRITE<6,77)  XI , X2 , X4 , X5 , FI , F2 , F4 , F5 , PI , P2 , P3 , P4 , P5 , P6 . 

C  1  P7,P8.P9,P10 

77  FORMAT <5F20. 8) 

50  DO  1400  II=1,N 

100  CALL  LRNDPC (DSEED,  U,  1) 

U  =  U*P10 

IF  <U. GT. P4)  GO  TO  500 
IF  <U. GT. PI)  GO  TO  200 
X  =  X2  +  U/F2 
GO  TO  1400 
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200  IF  (U.GT.P2;  GO  TO  300 
X  =  X3  +  (U-P1)/F4 
GO  TO  1400 

300  IF  (U.GT.P3)  GO  TO  400 
X  =  XI  +  (U-P2)/F1 
GO  TO  1400 

400  X  =  X4  +  (U-P3)/F5 
GO  TO  1400 

500  CALL  LRNDPCCDSEED,  V,  1) 

IF  (U. GT. P5)  GO  TO  600 
X  =  X2  +  (X3-X2)*V 

IF  < <U-P4)/(P5-P4)  .LE.  W)  GO  TO  1400 
V=F2  +  <U-P4)  /  <X3-X2) 

GO  TO  1300 

600  IF(U.GT.P6)  GO  TO  700 
X=X3+<X4-X3)*V 

IF  < (P6-U)/(P6-P5)  .GE.  W)  GO  TO  1400 

V  =  F4  +  (U-P5)/ (X4-X3) 

GO  TO  1300 

700  IF  (U  .GT.  PS)  GO  TO  900 
CALL  LRNDPC<DSEED,  V2,  1) 

IF  <V2  .GT.  V)  V=V2 
IF  <U. GT. P7)  GO  TO  800 
X=X1+  (X2-X1)J*:W 
V=Fl+2. *V*(U-P6)/(X2-X1) 

IF  <V.LE.F2l:V)  GO  TO  1400 
GO  TO  1300 
800  X=X5  -  W1:(X5-X4) 

V  =  F5+2. <U-P7)/(X5-X4) 

IF  <V.LE.F4*V)  GO  TO  1400 
GO  TO  1300 

900  IF  <U.GT.P9)  GO  TO  1000 
U=  (P9-U)/<P9-P8) 

X=X1-AL0G<U)/XLL 
IF  <X. LE. 0)  GO  TO  100 

IF  (V.LT.  (XLL* (Xl-X)+1. )/U)  GO  TO  1400 

V  =  W*F1*U 
GO  TO  1300 

3000  U  =  (P10-U)/(P10-P9) 

X=X5-ALOG(U)/XLR 

IF  (V.LT.  (XLR*(X5-X)  +  l.)/U)  GO  TO  1400 

V  =  V1:F5»U 

1300  IF  (ALOGCV)  .GT.  X3;*:ALOG<X)  +  R3X  -  X)  GO  TO  100 
1400  G(II)=X 
RETURN 

C  -  GAMMA  GENERATOR  A<=1.  ALGORITHM  GS. 

1500  CONTINUE 

IF  (A  .EQ.  ASAVE)  GO  TO  1550 

ASAVE  =  A 

AINV=1/A 

B=A*. 36787944  +  1. 

1550  CONTINUE 
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IF  <A  .EQ.  1.0)  GO  TO  1800 
DO  1700  1  =  1, N 
1560  CONTINUE 

CALL  LRNDPC(DSEED, UN,2) 

EXPL=-ALOG(UN(l) ) 

P=B*UN(2) 

IF(P  .GT.  1.)  GO  TO  1600 
X=P**AINV 

IFCX  .GT.  EXPL)  GO  TO  1560 
G(  I  )=X 
GO  TO  1700 
1600  CONTINUE 

X=-ALOG( <B-P)*AINV) 

IFC  <l-A)*ALOG (X)  .GT.  EXPL)  GO  TO  1560 
G(I)=X 

1700  CONTINUE 
GOTO  1900 
1800  CONTINUE 

CALL  LRNDPCCDSEED,  G,  N) 

DO  1810  1  =  1, N 
G(I)  =  -ALOGCGd)) 

1810  CONTINUE 
1900  CONTINUE 
RETURN 
END 
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